71 4538 


DASA  2475 


MELT- DOMINATED  IMPULSE  EXPERIMENTS 
AND  CALCULATIONS 


FINAL  REPORT 


September  1970 


J.  Reaugh,  A.  Lutze,  and  G.  Yonas 


f ; . 


Prepared  for 

DEFENSE  ATOMIC  SUPPORT  AGENCY 
Washington,  D.C.  20305 

Under  Contract  DASA-01-69-C-0095 
NWER  Subtask  4C321 


o 
.  c 


Preparing  Agency 

Physics  International  Company 
2700  Merced  Street 
San  Leandro,  California  94577 


Rnpre-dufCil  bY  _  .  , 

national  technical 

INFORMATION  SERVICc 

Vi  22I5I 


D  D  c 


/  Ainu  •  f  i  l 

l  Effl)  jli 


EE GDlTi 


'  +  ’*  **  HK.  J 

.  lot  —ItCM  Oid  u».  j 

la  uailmltad.  I 


DASA  2475 


PIFR-144 


MELT- DOMINATED  IMPULSE  EXPERIMENTS 
AND  CALCULATIONS 


FINAL  REPORT 


September  1970 


J.  Reaugh,  A.  Lutze,  and  G.  Yonas 


Prepared  for 

DEFENSE  ATOMIC  SUPPORT  AGENCY 
Washington,  D.C.  20305 

Under  Contract  DASA-01-69-C-0095 
NWER  Subtask  4C321 


Prepared  by 


Physics  International  Company 
2700  Merced  Street 
San  Leandro,  California  94577 


ABSTRACT 


A  combined  experimental  and  analytical  study  of  impulse  genera¬ 
tion  in  aluminum  has  beer,  carried  out  using  intense  electron  beams 
as  the  source  of  high  energy  density  loading.  Both  analytical  models 
and  hydrodynamic  codes  were  used  in  the  modeling  of  the  observed 
material  response.  In  addition  this  report  presents  a  detailed 
survey  of  current  calculational  techniques  and  an  extensive  dis¬ 
cussion  of  the  experimental  methods  required  for  this  work.  For 
the  electron  beam  fluence  range  employed  here,  giving  peak  doses 
less  than  2000  cal/g,  a  model  of  liquid  spall  was  found  to  represent 
the  data  well.  This  model  was  simplified  to  an  analytical  expres¬ 
sion  which  depends  exponentially  on  deposition  time  and  is  suitable 
for  predicting  melt-dominated  impulse  in  a  wide  range  of  materials . 

In  addition,  hydrodynamic  code  calculations  have  been  carried  out 
and  comparisons  are  made  with  measured  transmitted  stress  histories 
and  liquid  ejecta  velocities. 
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Manager,  and  Messrs.  Andrew  Lutze  and  John  Reaugh,  Program  Physi¬ 
cists  . 
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SECTION  I 


INTRODUCTION 


For  several  years  intense  electron  beams  have  suggested 

as  a  suitable  laboratory  too]  for  investigating  the  high  energy  density 
response  of  materials.  Recent  developments  in  beam  handling  and 
diagnostic  techniques,  along  with  improvements  in  laboratory 
methods  of  diagnosing  material  response,  have  permitted  workers 
in  this  field  to  collect  data  that  can  be  of  great  value  in  deter¬ 
mining  material  properties  and  in  evaluating  the  hydrodynamic 
codes  or  analytical  methods  that  use  these  material  parameters. 

The  essence  of  this  technique  is  to  use  a  well-characterized 
electron  deposition  profile  with  a  sufficiently  large  area  of 
irradiation  that  one  dimensional  stress-wave  propagation  can  be 
assured  in  the  solid  material  and  one-dimensional  flow  in  any 
material  ejected  from  the  front  surface.  In  addition,  one  is 
interested  in  achieving  high  enough  energy  densities  such  that 
material  phase  changes  can  occur.  To  accomplish  this,  one  must 
work  with  relatively  high-current  electron  beams  using  low- 
impedance  pulsers.  With  the  advent  of  low-impedance  transmission 
lines  and  diode  structures  between  1965  and  1967,  it  became  possible 
to  generate  the  required  beams .  Beam  physics  research  has  been 
continuing  since  that  time  and  beam  control  and  diagnostic  methods 
have  been  constantly  evolving  (References  1-3).  In  addition, 
pulsers  capable  of  generating  increasingly  higher  currents  have 
been  developed  and  the  present  state  of  development  should  be  viewed 
as  an  intermediate  one  (Reference  4)  . 

The  goal  of  this  work  was  to  investigate  the  response  of  a 
relatively  simple  material  to  high  energy  density  loading  in  order 
to  provide  a  firm  empirical  basis  for  analytical  modeling.  The 
state  of  development  of  the  electron  beams  at  the  time  of  this  work 
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Wet  -such  that  peaK  doses  less  than  2000  cal/g  could  be  achieved 
O'.,  r  ureas  of  3  cn"  using  the  Physics  International  738  Pulserad. 

Such  a  beam  was  therefore  used  to  study  the  process  of  impulse 
generation  resulting  from  blow-off  and  the  material  chosen  was 
aluminum.  Aluminum  itself  was  used  as  a  sample  materia]  for 
several  reasons:  (1)  it  is  a  relatively  easy-to-obtain,  inexpen¬ 
sive,  representative  metal,  (2)  it  has  no  major  solid-solid  phase 
transitions,  nor  is  it  particularly  rate  sensitive,  so  that  the 
calculation  of  stress  wave  propagation  is  straightforward.  The 
results,  however,  are  thought  to  apply  to  all  homogeneous  materials 
that  undergo  a  solid-liquid  transition  in  the  energy-density  range 
of  interest. 

Since  the  completion  of  the  experimental  portion  of  the  program, 
further  techniques  have  been  developed  using  the  738  Pulserad,  which 
allow  one  to  achieve  doses  up  co  6000  to  7000  cal/g  uniformly  over 
areas  <->f  roughly  0.5  cm  (References  5  and  6)  and  it  is  expected 
that  the  extension  of  the  work  discussed  here  to  this  higher  dose 
range,  will  be  relatively  straightforward.. 

This  report,  then,  presents  an  in-depth  approach  to  one  of  the 
fundamental  problems  characteristic  of  impulse  generation;  namely, 
that  dominated  by  removal  of  material  in  a  liquid  phase.  In  ad¬ 
dition,  the  report  presents  a  rather  general  discussion  of  the 
analytical  and  hydrodynamic  computer  code  modeling  of  impulse 
generation  over  the  full  rangy  of  phenom  nology  from  solid  spall 
to  full  vaporization.  Finally,  suggestions  are  presented  as  to 
how  this  method  can  be  extended  to  cover  the  entire  scope  of 
impulse-generation  effects. 
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SECTION  II 


ANALYSIS  OF  IMPULSE  GENERATION 


A.  INTRODUCTION 

* 

This  report  is  concerned  with  the  generation  of  a  net  impulse 

in  a  material  by  sudden  deposition  of  energy  through  the  use  of  a 

high-intensity  electron  beam.  The  process  of  energy  deposition  by 

electrons  slowing  down  in  a  material  in  itself  imparts  essentially 

no  momentum.  (At  the  nominal  electron  beam  parameters  used  in  the 

2 

experiments,  50  cal/cm  of  0.25  MeV  electrons,  the  total  beam 
momentum  is  less  than  0.01%  of  the  typically  measured  values.) 

Since  momentum  is  conserved  in  the  motion  of  deformable  bodies, 
the  imparting  of  a  net  impulse  requires  ejecting  material  from  the 
front  (irradiated)  surface.  The  content  of  this  report  is*  essen¬ 
tially  a  study  of  the  way  in  which  material  is  ejected. 

Before  presenting  a  detailed  discussion  of  calculational 
techniques  for  predicting  impulse,  it  will  be  useful  to  give  a 
qualitative  picture  of  the  processes  involved  in  removing  material 
from  a  suddenly  heated  region.  As  a  benchmark  for  this  discussion, 
we  will  consider  the  case  where  no  net  impulse  is  generated — a  solid 
that  stays  solid. 

Suppose  a  finite  region  inside  a  material  is  heated  instan¬ 
taneously.  There  is  then  a  sudden  localized  temperature  rise, 

T  -  T  .  Since  tne  material  is  compressible,  it  takes  time  to  expand 
to  the  relaxed  volume  VQ(1  +  3a (T  -  T  ) )  where  a  is  the  linear 
coefficient  of  thermal  expansion.  The  material  is  effectively 

* 

As  used  throughout  this  report,  impulse  is  defined  as  momentum 
per  unit  area. 


.  omprc"  ?ed  by  a  thermal  strain  in  the  amount  3a(T  -  T  ).  This, 
t,  v>n  ceases  a  stress  P  =  K(3u)  (T  -  T  )  ,  where  K  is  the  bulk  modulus. 
If  e  ie  the  specific  deposited  energy,  the  formula 


P 


KIM-  e 

C  e 


holds,  where  Cv  is  the  specific  heat  and  y  is  the  Gruneisen  ratio. 
If  we  assume  that  the  material  is  elastic  and  that  conditions  of 
one-dimensional  strain  hold,  the  stress  pulse  will  separate  into 
two  elastic  waves  moving  in  opposite  directions  with  half  the 
original  amplitude.  If  the  heated  region  starts  at  a  free  surface, 
one  (compressive)  wave  moves  into  the  cold  material.  The  wave 
moving  initially  toward  the  heated  surface  is  reflected  with  op¬ 
posite  sign  and  moves  into  the  solid.  In  this  way,  then,  a 
suddenly  heated  region  near  the  front  surface  generates  a  compres¬ 
sive  and  tensile  stress  pulse  moving  into  the  material.  A  point 
inside  the  material  will  feel  first  the  compressive  (direct)  part, 
then  the  tensile  (reflected)  part.  In  this  way  there  is  no  net 
impulse: 


CO 


O 


As  was  stated  earlier,  in  order  to  generate  a  net  impulse  by 
sudden  heating,  material  must  be  ejected  at  the  front  (irradiated) 
surface.  This  ejection  is  basically  one  of  two  types — the  motion 
of  spalled  or  vaporized  material  away  from  the  solid.  Much  of  the 
difficulty  with  describing  impulse  generation  quantitatively  is 
associated  with  the  fact  that  the  two  types  of  ejection  are  not 

J  i~  *  s*  4-  •->  vn  t*vn1  u  nnltf  annrnnri  af  o  »  4*  fho  oyfronribC  nf  f’Am- 

U.,0  UAAIWW./  UtiW*  WA.  C''*'  WW-  *-  —**  — ... - —  —  —  - - - 

perature  rise. 
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°.pall,  the  low- temperature  phenomenon,  is  simply  tensile 
failure.  The  (reflected)  tensile  stress  wave  exceeds  the  cohesive 
strength .  If  this  occurs  at  only  one  depth  in  the  material,  all 
the  material  moving  away  from  the  solid  remainder  is  ejected;  if 
it  occurs  over  a  region,  that  region  solid  is  pulverized.  For  the 
case  at  hand,  the  spalled  material  has  a  uniform  velocity  away  from 
the  solid,  since  that  is  the  direction  of  material  motion  in  the 
reflected  wave.  The  fact  of  spall  implies  that  the  stress  wave  in 
the  remaining  solid  will  no  longer  be  symmetric — the  tensile  part 
will  no  longer  exceed  the  cohesive  strength — and  a  net  impulse  is 
generated. 

On  the  other  hand,  very  large  temperature  rises  will  produce 
vaporization.  The  vapor  state,  unlike  the  solid,  has  the  property 
that  it  is  never  in  tension.  The  fact  that  it  is  never  tensile, 
however,  implies  that  there  is  also  no  restoring  force  once  it 
begins  to  expand.  As  a  consequence,  the  vapor  moves  away  from  the 
solid,  and  the  stress  pulse  in  the  solid  is  a  compressive  pulse 
that  decays  slowly  as  the  pressure  of  the  vapor  on  the  solid  sur¬ 
face  decays. 

Neither  of  the  above  pictures  of  behavior,  which  may  be 
considered  a  mechanical  description,  are  totally  accurate.  This 
is  seen  to  be  the  case,  particularly  when  one  considers  a  phase 
diagram  drawn  from  considerations  of  equilibrium  thermodynamics. 

The  paths  to  be  considered  on  the  equation-of-state  surface  are 
the  initial  pressure  (and  temperature)  rise  at  constant  volume, 
and  the  adiabatic  release  paths . 

In  Figure  1,  in  particular,  two  release  adiabats  are  shown 
on  such  a  phase  diagram.  The  heavy  line  represents  the  path 
taken  by  the  spall  model,  used  in  this  case  to  describe  liquids. 


FIGURE  1.  SCHEMATIC  LOADING  PATH  AND  RELEASE  ADIABATS 


The  solid-liquid  phase  transition  is  ignored,  and  the  liquid  is 
assumed  to  release  to  just  below  zero  pressure,  the  liquid  cohesive 
strength,  followed  by  spalling  of  the  liquid  so  that  zero  pressure 
is  maintained.  The  actual  path  taken  is  probably  more  like  the 
dashed  line.  The  liquid  release  adiabat  is  followed  down  to  some 
positive  pressure,  the  metastable  limit.  It  then  jumps  back  to 
some  equilibrium  vapor  pressure,  and  further  expansion  occurs  in¬ 
side  the  liquid-vapor  dome.  Thu  liquid  spall  model,  however,  is 
found  to  describe  the  behavior  of  liquid  spall  quite  well.  This 
can  only  be  the  case  when  the  equilibrium  vapor  pressure  is  quite 
small  with  respect  to  the  initial  pressures,  so  that  ignoring  the 
vapor  pressure  affects  the  solution  very  little.  For  solid  spall, 
of  course,  the  equilibrium  vapor  pressure  is  less  than  the  triple 
line  pressure  (for  materials  that  exhibit  a  liquid  state  it  is 
less  than  1  atm)  so  that  the  pressure  can  be  considered  zero. 

Vapor  expansion  adiabats,  too,  will  cross  into  the  liquid- 
vapor  coexistence  region  (Reference  7) .  If  the  initial  deposited 
energy  density  is  high  enough,  however,  it  will  occur  at  sufficiently 
low  pressures  so  as  not  to  materially  affect  the  impulse  generated. 

There  are  several  methods  of  calculating  the  impulse  generated 
by  rapid  energy  deposition  in  a  slab  of  material  at  ambient  con¬ 
ditions.  These  methods  may  be  grouped  into  two  categories,  that 
use  respectively,  "analytical  expressions"  and  "hydro-codes."  The 
analytical  expressions  require  the  evaluation  of  a  one-dimensional 
integral,  which  is  the  total  impulse.  The  hydro-codes  are  of  a 
family  of  computer  programs  evolved  from  programs,  developed  at 
the  University  of  California's  Lawrence  Radiation  Laboratory,  that 
require  an  equation-of-state  of  the  form  P(V,E).  These  programs  are 
used  to  calculate  the  detailed  motion  (typically  in  one  dimension) 
of  the  expanding  material.  The  total  impulse  is  then  the  asymptotic 
value  of  this  momentum  at  long  time  intervals. 
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B. 


ANALYTICAL  EXPRESSIONS 


Two  expressions  are  in  use  that  give  impulse  as  an  integral. 
They  are  the  Whitener  formula  and  the  BBAY  formula.  In  their  most 
general  form,  they  are  given  in  Equations  (1)  and  (2) ,  respectively. 


a.  Whitener  Formula 

The  Whitener  formula  may  be  derived  from  consideration  of 
an  incremental  mass  inside  the  expanding  material.  By  con¬ 
sidering  the  final  state  energy  to  be  partitioned  into  the  kinetic 
energy  of  the  hydrodynamic  motion,  and  the  final  state  internal 
energy  density, 

h  Am  (V2)  +  Ef  Am  =  Am  (3) 

where  E^  is  the  initial  energy  density.  Obtaining  the  expression 
for  momentum,  by  assuming  unidirectional  velocity,  and  summing  the 
contributions  gives 

I  =  ^  V  dm  =  /2  y*o/TB7“="EJ)r  dm  (4) 
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b .  EBAY  Formula 


The  BBAY  formula  (Reference  8 )  is  a  closed  form  solution  to 
the  one-dimensional  equations  of  motion  and  conservation.  The 
hydrodynamic  equations  are 


3x 
3  c 


V 


3x 

3t 


U 


3U  3£ 

3t  "  3o 


3e 

3t 


D  (a,  t) 


(5) 


where  D(c,t)  is  the  heating  rate,  x  is  the  Eulerian  (space)  co- 

X 

ordinate,  x  is  the  Lagrange  (material)  coordinate, o ^ po(x')dx' 
is  the  areal  density  from  front  surface,  V  is  volume,  U  is  material 
(particle)  velocity,  and  e  is  the  internal  energy  density.  The 
fourth  expression  contains  the  hydrodynamic  assumption  of  small 
distortional  energy  and  no  viscosity. 


It  is  assumed  that  the  deposition  occurs  at  a  constant  rate, 
and  that  the  final  internal  energy  is  subtracted  out,  so  that 

D(a;t)  =  |  max  [E^c)  -  Ef  (a)  ,0]/pQ (a)  tx ,  0  £  t  <  t1 
'  ®  t^  <  t 

where  E^(a)  is  the  total  energy  per  unit  volume  to  be  deposited  at 

w  •  xiic  ^uuikvxujr  #iuo  uic  ujr  uaau 


W  “  W 
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Defining  a  new  dependent  variable  Z  by 

ra 

Z  =  j  x(a’,t)d0* 

the  expressions  of  Equation  (5)  become 

X  =  Z ' 

V  =  Z '  ' 

u  =  z  * 

P  =  -z 

e  *  D(0;t)  +  Z  Z*  '  (6) 

If  the  equation  of  state  of  the  vapor  is  given  by 

e  =  pV/(y~l) 

the  y-law  gas  as sumption , 

e  =  (pV  +  pV)/(y-l). 

Thus  the  e  equation  f  (6)  becomes 

D(0;t)  +  Z  Z"  =  [-  Z  Z"  -  Z  Z "]/ (y-1)  , 

or 

Z  Z"  +  yZ  Z"  =  -  (y-1)  D  (o,t)  . 

Assuming  that  Z  is  separable  t 

Z(o,t)  =  Y (0)  f(t)  ,  and 
[f  f’+  yf  f ] YY"  =  -  (y-1)  D(0,*t)  . 
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Since  D  is  a  function  of  o  only, 


and 


f  f  +  Yf  f 


r  ,  o  ^  t  <  t, 

s  1 

0 ,  s:  t 


(7) 


Y  Y"  =  -  (Y-l)  D(o,  t)/r  =  -  (Y-l)max  [E.  (o)  -  E- (o)  ,  0]  /  (I’  p  (a)  t.)  (8) 

i  L  SOI 

where  rg  is  the  separation  constant. 

Equation  (7)  is  solved  by  f  =  tm,  0  £  t  <  t, ,  which  becomes 


t2m  3  [m(m-l) (m-2)  +  Ym2  (m-1)] 


T  implying 

b 


"  -  I  '  rs  "  I  (3^-1> 


2m  -  3  =  0  so 

t  £  tv  f"(Y_1)  (fY  f)  =  0,  or 


For 


f*J  =  const  =  t,3^  3  (^V1/2  =  3  tl{3Y-1}/2 

Thus  for  t  >  t1  and  using  f(t1)  =  t^2 ,  one  obtains  the  solution 

■3t.  3t,  (3Y"D/2  ,  .  .1  k 


rjti  jti 

f  [23  "  2TF 


IT 


:-(Y- 


fl)l  ^ 

/ 


where 


-  3  l  s (v-i ) 

P  “  L2  +  Y3!^  “  ~4T 

1  s 


Equation  (8)  becomes 


Y  Y"  =  -  - 


J  p- ('oTtT  "lax  [E1(0)  "  E2(o,'0]'  0  s  °v 
o  1 
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with  the  boundary  conditions 

P(0,t)  =  -Y (0) f  (t)  =  0 

U(ov)  =  Y'  (av>f  (t)  =  0 

That  is,  the  left  bour.aary  is  a  free  surface,  and  the  solid-vapor 
interface  is  (nearly)  fixed. 


The  impulse  is  defined  as 

*8  =  X  PlVt)dt 

/•»  .. 

=  -  /  z  dt 

Jo 

=  Z(av,0)  -  z (av,«) 

=  [f  (0)  -  f(«)]Y(av) 


Since 


ra 

f(0)Y(o)  =  Z(0,a)  =  /  U(o',0)da'  and  since 


U(a,0)  =  0  and  f(0)Y(o)  j-  0  implying  Y  (a)  ^  0, 

f  (0)  =  0. 

Since  the  analytical  form  of  f  is  a  monotonically  increasing 
function,  (without  boundary) 

f  (°°)  = 
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Defining 


F (a)  =  t  - - -  max  [E,  (a)  -  E_(c)  ,  0|  , 

po(0)tl 

Equation  (8)  with  U(c?v)  =  0  B.C.  implying  Y'(ov)  =  0  becomes 


Y'(o>  -  -  /  T&Td' 

•Srf 

°V 

and  P  (0 ,  t)  =  0  B.C.  implying  Y { 0 )  =0  gives 

a  r  °v 

Y  la)  =  f  do"  I  f 


a 

A  do" 

r  /V  F(o,)  do- 
/  Y(O-) 

r 

0 

*/a" 

For  o  ^  a  one  can  interchange  the  order  of  integration,  since  the 
integral  is  not  a  function  of  o’’  , 

a.. 


Y  <Ov)  - 


fv  g'Ffo': 
J  Y(0r 


Equation  (9)  can  be  solved  numerically  by  an  iteration  technique, 

but  the  standard  BBAY  equation  states  that  the  average  value  of  Y 

2 

over  the  interval  0  to  a is  given  by  Y(o  )/a  .  Thus  substituting 
the  average  for  Y(o'), 


?  2  f  v 

[  Y  (cfv)  3  «  o  /  O '  F  { a ' )  do '  , 

o 

=  “  V2F"  ^/t  }/  da'  Po(0')  maX  ^Ei(°,)“Ef  (o')  ,Cip 


or  simplifying 


r  o 

yb  -  “[2  /  [Ei!£,,)  ■  E£(0’)!  d“'  ^  <11 

_  2 
It  is  then  shown  that  for  reasonable  forms  of  E^(c) ,  I  ^  a  <  2; 


hence  a  is  chosen  to  be  1.2. 
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The  form  of  the  BBAY  formula  (although  not  the  coefficient) 
may  be  derived  by  considering  a  snowplow  model.  At  a  point  a 
distance  x  into  the  expanded  material,,  which  corresponds  to  m  mass 
units,  the  velocity  of  the  Am  mass  is  lower  than  in  the  Whitener 
form,  since  it  pushes  on  the  mass  ahead  of  it.  Thus  the  kinetic 
energy  is  h  mv2  =  (E^  -  Ef)Am,  and  the  square  of  the  momentum 
contribution  at  a  point  is 

(mv)2  =  2m(E^  -  Ef)Am 

Thus  the  total  impulse  squared  is 

x 

I2  =  f  2m  (E.  -  E  )  dm 
and  the  impulse  is  just 


which,  except  for  the  coefficient  a  is  the  BBAY  formula. 

2 •  Use  of  Analytical  Expressions 

Both  the  Whitener  and  BBAY  formulas  contain  an  adjustable 
parameter  Ep(x)  that  corresponds  to  the  internal  energy  of  the 
final  state.  There  are,  in  general,  two  ways  in  which  this  may 
be  obtained. 

a .  Final  Energy  as  Polynomial  of  Deposited  Energy 

The  first  way,  which  has  had  some  success,  is  to  plot  experi¬ 
mental  impulse  data  and  a  family  of  lines  of  calculated  impulse 
with  Et-  the  parameter.  This  produces  a  plot  such  as  shown  in  Figure  2. 
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FIGURE  2.  TYPICAL  CALCULATED  IMPULSE  WITH  Ef  PARAMETER 
AMD  IMPULSE  DATA 


Using  such  data,  one  can  obtain  a  least-squares  fit  to  the 
final  energy  density  as  a  function  of  the  initial  (deposited)  energy 
density,  assuming  a  polynomial  form.  For  each  data  point,  the  depth  of 
material  removal  is  measured  (or  estimated  by  the  temperature  at 
which  the  tensile  strength  goes  to  zero) .  Assuming  that  this  is 
constant,  or  taking  an  average  value,  one  has  the  minimum  final 
energy  Eq.  Then  for  each  data  point,  one  has 


=  K  j 
Jo 


x(E*) 


X(E.  (x)  -  E*)dx  = 


■/ 

Jo 


X(EQ)  r 


X 


(Ei(x) 


Ef  (e.  U) ) 


dx 


where  E*  is  the  apparent  final  energy,  obtained  by  interpolating 
impulse  between  lines  of  constant  Ep.  Rearranging, 


x(Bj 


f  °  xEi(x 

•4(e*) 


)dx  -  e*  LxJ.|*Ll 


x(E_) 


:h 


f  °xVEi 

Jet 


.  (x) ) dx 


For  the  k'"  data  point,  call  the  left-hand  side  (LHS)  -D^. 
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Assuming 


-  £  “i  Ei  > 


where 


a  =  E 
o  o 


and  defining  for  the  1c  data  point 


x(E  ) 


r  O 

/  x[Ek(x)]n  dx  ,  n 

J  n 


0,  m  , 


i.e.,  n+1  quantities 


=  -  V*  a.  J. 
Z-f  3  3k 


Using  the  standard  definition  of  a  least-squares  fit,  the  system 
of  m  equations  (m  <  N,  the  number  of  data  points)  for  the  eu  is 


(M±j)  («j)  =  (PjJ  /  i, j  -  ltm  , 


where  the  elements  of  column  vector  3.  are  (a  =  E  ,  3  ,  =  0.5  x^(E  )) 

1  '0  o'  ok  o 


‘i  =  E  (Dk  -  Eo  Jok>  Jik  ' 


and  the  elements  of  the  symmetric  m  by  m  matrix  Mi;.  are 


m..  =  r, 

1]  Jt-J  X 


ik  Jjk 
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Thus  the  solution,  (symbolically) , 


is  given  by 


a , 


(Mi.) 


-1 


(V 


This  solution,  of  course,  does  not  require  that  the  data 
points  have  the  same  deposition  profile,  so  that  data  from  widely 
varying  sources  may  be  used  together,  provided  only  that  E*  and 
x(E*)  can  be  obtained,.  For  this  least-squares  fit  to  be  useful 
over  large  ranges  of  energy,  large  amounts  of  reliable  data  or 
considerable  faith  in  extrapolation  is  required,  for  the  value 
of  impulse  is  certainly  related  to  the  value  of  E^. 

b.  Freezing  (McCT)  Model  for  Final  Energy 

The  alternative  method  for  obtaining  Ef,  and  the  way  required 
when  data  are  lacking  or  uncertain,  is  to  construct  a  model  for 
the  final  state  of  the  partially  condensed  gas,  and  derive  a 
functional  form  for  Ef  (E. ) .  The  most  successful  of  the  models 
has  been  the  McCloskey-Thompson  model.  This  model  is  discussed 
in  some  detail  below. 


It  is  known  (Reference  7)  that  at  some  point  in  the  isentropic 
expansion  of  gas,  the  composition  freezes.  The  distance  between 
the  condensing  droplets  becomes  large,  requiring  a  longer  interval 
to  reach  thermodynamic  equilibrium.  When  this  time  interval 
becomes  long  with  respect  to  the  time  associated  with  hydrodynamic 
motion,  condensation  stops,  and  the  composition  remains  fixed. 
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The.  McCloskey-Thompson  (McCT)  model  assumes  that  the  composition 
stays  fixed  at  the  intersection  of  the  isentrope  with  the  triple  line. 
This,  of  course,  is  a  useful  assumption,  since  pressure  and  tempera¬ 
ture  are  constant  on  the  triple  line  as  the  specific  volume  changes. 
Their  derivation  then  proceeds  as  follows: 


Consider  a  solid  heated  instantaneously.  The  initial  energy 
and  entropy  are  given  by 


E.i  "  E«.  <VTt>  +  cv>.  (Ti  -  V 


si  “  s*  <W  +  cvl  ln 


ft) 


for  CvJl  a  constant.  Denoting  Eq  and  SQ  as  the  constant  (base) 
energy  and  entropy,  and  evaluating  them  at  the  normal  melting 
temperature  (which  is,  in  general,  near  the  triple  temperature) 
one  can  solve  for  (E^)  by  eliminating  in  Equation  (11) 
and  obtain 


Si.  ”  So  +  cvl 


/  B.-E  \ 

^  + 


If  the  final  state  is  on  the  triple  line, 


sf  =  xSo  +  (1_x,svo 
E(  =  XE0  +  (l-x)Evo  . 


(13! 
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Eliminating  x,  the  liquid  fraction,  and  solving  for  the  final 
entropy, 


since 


■  so  +.^r  <Ef  -  V 

VO  o 


AS 

AE 


*=£L 

1-9 


AS  _  1_ 

AH  ~  T 

o 


that  is,  the  change  in  enthalpy  for  a  constant  pressure  phase 
transition  is  nearly  the  change  in  internal  energy,  the  final 
entropy  is  given  by 


S 


f 


V  * 


Equating  the  initial  and  final  entropies  for  the  isentropic 
expansion, 


Ef  =  Eo+Cv*T 


Approximating  C.^  Tq  to  E^,  the  final  form  is  given  by 


('  E.  -  E  \ 

1  + 


f  o 


1  +  In 


ft)]- 


The  maximum  initial  energy  for  which  this  formula  is  valid  is 
given  by 


vo 


=  E 


(1  + 
°\ 


Jin 


E.  „  \ 

•j  — m  v  (?  \ 


or 


(14) 
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(15) 


E .  ,,  =  E  exp 

1  maxi  o  ^ 


E  -  E 
vo  o 

E 


and  the  minimum  energy  Eq. 


e.  Non-Freezing  Model  for  Final  Energy 

One  may  reinterpret  the  McCT  model  and  remove  the  essentially 
arbitrary  choice  of  freezing  the  composition  at  the  triple  tem¬ 
perature  and  pressure.  For  materials  that  are  in  a  liquid  state 
at  ambient  pressure,  the  pressure  of  the  triple  point  is  less  than 
atmospheric — water,  for  example,  and  a  triple  point  pressure  of 
5  Torr.  Most  impulse  experiments  are  performed  at  1  atm  or  less. 

It  is  improbable  that  impulse  measurements  made  at  1  atm  or  less 

* 

will  change  with  ambient  pressure.  Thus  we  may  replace  the 
assumption  that  the  composition  freezes  by  the  assumption  that 
the  material  expands  to  ambient  pressure,  which  is  taken  to  be 
the  triple  point  pressure. 


Since  material  vapors  generally  follow  the  ideal  gas  law  at 
pressures  less  than  1  atm,  one  can  easily  calculate  the  final 
entropy  and  internal  energy . 


For 


E.  >  E. 
l  l  max 


Ef  =  Evo  +  Cvg  (Tf  -  V 

/T-\  /V 

Sf  -  Svo  +  Cvg  *n  \tl+R  *nl  V 

\  o  \  o, 


(16a) 

(16b) 


where  by  applying  the  assumption 


Pr  =  P 
f  o 


RT 


v, 


Considering  that  most  of  the  impulse  is  generated  in  1  nsec  or  so, 
1  bar-usec  corresponds  to  1  tap. 
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Setting  Evq  -  C  Tq  =  Egub  and  solving  Equation  (16a)  for  Vf, 


V, 


R (E-j  -  E  ,  ) 
f  sub 

C  P 
vg  o 


Then  Equation  (16b)  becomes 


S.  =  S  +  C  in 
f  vo  vg 


E,  -  E  . 
f  sub 

C  T 
vg  o 


+  R  in 


E-  -  E  . 
r  sub 

C  T 
vg  o 


(17) 


Equating  Equations  (17)  and  (12) ,  replacing  C  .  T  with  E  ,  and 

V  jC  O  O 

C  T  with  E  .  -  E  , 

vg  o  vo  sub 


in 


( 


Ef  Esub  \ 
Evo  Esub  / 


vi 

R3 


or 


Et  -  Esub  +  (E 


-  E 


vo  sub 


ft)- 

•*)(v  Exp 


s  -  s 

vo  o 


R3 


('«) 


where 


»  ■  ■ 


If  one  makes  the  additional  replacements  that  Egub  -  L, 

AS  =  L/T  ,  where  L  is  the  change  in  enthalpy  between  liquid  and 
vapor  at  tho  melt  temperature. 


E 


f 


Cvi/RB 


L  +  (E.„  -  L) 


VO 


ft)  expft^)- 


(18) 


By  construction,  Equation  (18)  has  the  property  that  Ef  (for 

T*  =  T  so  L  =  0  and  E.  =  E  )  =  E„  ,  so  that  it  joins  smoothly 
f  o  x  o  vo 

with  Equation  (15) . 
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Since  for  E-  >  E.  ,, ,  E,  is  an  increasing  function,  the  use 
i  x  max? '  f  J  ' 

of  a  constant  value  as  suggested  in  the  McCT  reference  will  give 
results  that  predict  a  larger  impulse  at  high  initial  energy 
densities . 

C .  HYDRODYNAMIC  CODES 


The  analytical  expressions  discussed  in  this  section  do  not 
deal  with  the  shape  of  the  stress  wave  propagated  into  the  remain¬ 
ing  solid  material.  If  this  stress  wave  shape  is  required,  one 
must  solve  the  hydrodynamic  equations  of  motion  and  conservation. 

In  addition,  as  was  seen  in  the  BBAY  derivation,  a  relation  between 
pressure,  volume,  and  internal  energy — an  equation  of  state — is 
required.  In  these  codes,  then,  it  is  the  equation  of  state  that 
determines  the  impulse  and  pressure  history,  so  that  a  study  of 
impulse  calculation  with  hydro-codes  is  a  study  of  the  various 
equations  of  state. 

1 .  One-Phase  Equations-of-State 

In  a  restricted  range  of  initial  energy  densities,  one-phase 
equations  have  been  used  successfully  in  predicting  the  impulse 
and,  in  fact,  the  transmitted  stress  profile  obtained  by  quartz 
gauge  measurements  (Reference  9) .  For  maximum  internal 
energy  densities  such  that  the  vapor  pressure  is  less  than  a  few 
bars,  equations  of  state  generated  from  Hugoniot  (shock  compression) 
data  with  two  modifications  are  used.  First,  for  elastic-plastic 
codes,  the  yield  stress  goes  to  zero  as  internal  energy  approaches 
the  enthalpy  of  a  melting  solid.  Second,  and  more  important,  the 
tensile  strength  goes  to  zero  as  the  energy  approaches  the  melt 
enthalpy.  If  this  tensile  strength  is  exceeded  at  a  point,  the 
pressure  is  set  to  zero  there,  and  the  material  is  presumed  to  have 
spalled  at  that  point.  In  this  way  a  net  momentum  is  imparted  to 
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the  solid  material;  the  transmitted  stress  profile  is  nearly  all 
compressive. 

If  appreciable  vapor  is  formed,  however,  impulse  calculations 
with  a  one-phase  equation  of  state  result  in  predictions  of  impulse 
lower  than  that  measured.  It  is  apparent,  then,  that  for  higher 
energy  densities,  a  description  of  the  vapor  phase  is  important. 

2 .  Two-Phase  Equations  of  State 

Two-phase  equations  of  state  in  current  use  incln.de  the 
Tillotson  equation  (Reference  10) /  its  modification  (Reference  n) 
and  the  PUFF  equation  (Reference  12) .  These  have  the  following  in 
common:  (1)  there  is  a  line  in  internal  energy-volume  space  that 

separates  the  two  equations  of  state  (these  lines  are  compared  in 
Figure  3) ,  and  (2)  the  two  equations-of-state  forms  are  usually 
continuous  at  the  boundary  and  have  continuous  first  derivatives. 
The  equations  are  given  in  Table  I. 

The  PUFF  and  modified  Tillotson  equations  have,  among  other 

parameters,  one  called  E  .  The  standard  Tillotson  equation  uses 

s 

E  =0.  At  large  expansions,  (n  small)  the  dominant  term  of  both 
s 

equations  is  of  the  form  with  P  equal  to  a  constant  (E-Es) ,  the 

y-law  gas  form.  Thus  to  fit  experimental  data  would,  in  general, 

require  the  use  of  Eg  as  an  adjustable  constant,  as  was  the  case 

in  the  BBAY  and  Wbitener  expressions.  The  PUFF  and  modified 

Tillotson  equations  give  reasonably  similar  results  for  impulse 

when  the  same  value  of  E  is  used.  This  is,  of  course,  in  agree- 

s 

ment  with  the  results  of  the  BBAY  derivation,  which  showed  that 
impulse  was  independent  of  the  value  of  y  in  the  y-law  gas 
equation  of  state.  They  are  also  both  in  reasonable  agreement 
with  the  BBAY  calculations  using  a  constant  Eg. 
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FIGURE  3  PHASE  DIVISIONS  OF  VARIOUS  EQUATIONS  OF  STATE 
(IN  ALUMINUM) 


TABLE  1 


Name 

PUFF 


TILLOTSON 


MODIFIED 

TILLOTSON 


TWO-PHASE  EQUATIONS  OF  STATE 


Solid  Vapor 
x 


Expression 


x 


X 


Cy+Dy2+Sy3+GnE 


x  n[  H+  (G-H)  rr*]  |e-Es[  1-exp  [g~  i  (l-i)]]j 

[a+b^§-  ^-1^  jEn+Ay+By2 
aEn+jbF,n^|-  ~+l^  +Ay  exp  B  (l~) 

ho  -  tf] 

h  ft  ?  *  f] 


x  aEn+ 


exp 


En  +  Ay  +  By 


x  a (E-E  ) p  +  G  exp  R, 

b 


G  =  aE  n  +  b 

5 


(f-Vi) 

\on  / 


-l 


En+(Ay+By  )exp(BR), 


If  the  solid  is  heated  so  that  the  entropy  is  greater  than 

the  critical  point  entropy,  the  distinction  between  phases  becomes 

quite  arbitrary;  properties  change  continuously  between  the  "phases." 

Thus  a  two-phase  equation  of  state  makes  good  physical  sense  for 

this  high  energy  density  regime  since  pressure  changes  smoothly 

along  an  isen trope.  If  i  choice  were  necessary  among  the  three 

equations  offered,  the  modified  Tillotson  equation  would  be  best, 

as  the  solid  Grur.eisen  coefficient,  T,  varies  from  the  ambient 

value  of  approximately  2.0  to  the  asymptotic  Thomas-Fermi  gas  value 

of  0.5  as  the  internal  energy  increases.  Whether  the  form  chosen 

for  interpolation  at  reference  density  0.5  +  (F  -  0.5)/(E/Eo  +  1) 

is  correct  is  not  relevant.  The  equation  is  incorrect  where  the 

* 

interpolation  form  is  important. 

The  inaccuracy  of  the  two-phase  equations  is  greatest  where 
the  mixed-phase  (liquid-vapor)  region  is  important.  The  calculation 
of  pressure  in  this  region  is  difficult,  particularly  if  there  is  a 
requirement  for  thermodynamic  consistency  and  agreement  with  the 
limited  equation  of  state  data  that  is  available. 


ic 

Another  equation  of  state  with  a  temperature-dependent  Gruneisen 
was  used  in  Reference  14  to  fit  Hugoniot  data  on  porous  media  at 
pressures  up  to  several  megabars.  It  does  have  the  slight  dis¬ 
advantage  of  being  implicit  in  energy  (temperature  is  the  dependent 
variable) . 
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D .  MONK  CODES 


MONK  codes  (Reference  13)  generate  a  two-phaso  equation  of 
state  by  an  approach  that  differs  from  other  two-phase  equations. 
Instead  of  selecting  eguation-of-state  forms  that  join  smoothly,  the 
MONK  code  uses  separate  forms  for  each  phase  and  joins  them  by 
equilibrium  thermodynamics. 


The  vapor  equation  is  one  suggested  in  Reference  16  as  being 
useful  for  moderately  high  pressures  (hundreds  of  kilobars)  and  at 
near  liquid  densities.  It  is  a  virial  expansion  cf  the  form 


P 


RT 

V 


where  B (T)  is  the  second  virial  coefficient  of  Lennard-Jones  (6-12) 
potential  and  C,  D,  and  E  are  constants.  The  coefficients  are 
given  as  dimensionless  ratios  of  the  co-volume.  The  temperature 
and  volume  scaling  parameters  are  related  to  the  two  constants  in 
the  6-12  potential,  which  are  the  depth  of  the  potential  well  and 
the  radius  of  the  minimum  potential.  They  may  in  turn  be  estimated 
from  the  critical  point  temperature  and  volume. 


The  solid  form  is  separated  into  two  parts  corresponding  to 
nuclear  and  electronic  contributions.  The  latter,  of  course,  are 


only  dominant  at  high  temperatures ,  The  nuclear  form  uses  the 
compressibility  and  its  pressure  derivative,  the  coefficient  of 


thermal  expansion,  and  the  specific  heat  at  constant  volume  to 


obtain  the  nuclear  contributions  to  pressure  and  internal  energy 
as  functions  of  temperature  and  volume.  This  solid  form  is  assumed 
to  hold  throughout  the  solid  and  liquid  phases. 
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It  is  desired,  then,  to  make  a  transition  between  the  forms  at 
high  mass  densities  and  temperatures  and  to  make  it  thermodynamically 


consistent.  For  phase  changes  at  constant  volume  and  temperature, 
the  thermodynamic  potential  is  the  Helmholtz  free  energy.  Since 
the  expression  for  the  Helmholtz  free  energy,  F  =  E  -  TS,  can  be 


written  for  both  phases,  equating  and  Fg  defines  a  line  Tsuk  • 
For  the  MONK  code,  then,  the  transition  between  the  two  phases  is 


a  line  in  the  T  -  V  plane,  but  a  region  in  the  E  -  V  plane.  In 
this  region,  the  pressure  is  calculated  by  a  simultaneous  solution 
of  the  three  equations 


n  E  +  (1  -  n  )  E  =E 
g  g  g  s 


\  vg  +  (1  '  V  vs  "  v 


P  (T  .  (V)  ,  V  )  =  P  (T  .  (V)  ,  V  )  =  P 

g  sub  g  3  sub  s 


As  with  the  other  two-phase  equations,  the  formulation  is 
only  appropriate  for  high  energy  densities,  where  the  release 
adiabats  remain  in  the  vapor  phase. 


For  adiabats  that  intersect  the  mixed  phase  region  (see 
Figure  4),  however,  it  is  necessary  to  construct  the  boundary 
of  that  region. 
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Solid-Liquid 


FIGURE  4.  TYPICAL  PHASE  DIAGRAM 


E.  MIXED-PHASE  REGIONS 

1 .  Construction  of  the  Mixed-Phase  Boundaries 

The  boundaries  of  the  mixed-phase  region  (liquid-vapor)  may  be 
established  by  the  following  arguments. 


For  any  N-component  mixture  we  may  write  the  change  in  Gibbs 
free  energy,  G,  as 

N 

dT  +  >  '  ^  dn±  +  (19) 


dG  ~  3Vt  n  dP  +  37p  n 

A'ni  /P'ni  i=l  A/P,T  ,r\^± 


E3G  n 
9n±> 


now 
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an, 


Mi ' 


the  chemical  potential,  and  for  equilibrium  between  two  phases  of 
the  same  element  at  constant  total  mass,  N  =  2,  p-j  =  p2,  dr^  =  -dn2/ 
so  that  Equation  (19)  reduces  to  dG  =  VdP  -  SdT.  However,  along 
a  line  where  T  is  equal  to  a  constant  in  the  mixed-phase  region, 

P  is  constant,  so  dG  =  0. 


Thus,  at  temperature  T1,  we  know  that  Gsoii^  V^) 

=  G  (T. ,  V_)  (the  subscripts  refer  to  Figure  5) .  The  gas  and 
93s  X  z 

solid  energies  may  be  placed  on  the  same  scale  by  assuming  that 
AH,  P^,  and  are  known  at  a  boiling  point.  Then  since 
AH  =  AE  +  A (PV) ,  and  for  standard  boiling  point  data,  P  is  constant, 
AHb  =  AEb  +  PfaAV. 
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Vg  is  found  from  the  expression  (T^,  V^)  =  P^;  Vg  is  found 
similarly . 

o 

If  E  can  be  represented  functionally  as  (T,V)  +  E 

AHb  =  VV  V  +  E°  -  Es  <V  V  +  pb  (vg  '  V  ' 

or 

E°  AH  -  Eg  <Tb,  Vg)  +  Es  (Tb,  Vs)  -  Pb  (Vg  -  Vs>. 

The  gas  and  solid  entropies  are  placed  on  the  same  scale  by  the 
equality  of  Gibbs  free  energy  at  the  boiling  point. 

G  =  F  +  PV  =  E  -  TS  +  PV 


again  write  the  gas  entropy  as  (T,  V)  +  S  . 

Eg  <V  V  +  E°  ‘  Tb  CSg  <V  V  +  S°]  +  Pb  Vg  ’ 
Es  (Tb'  V  -  Tb  [Ss  <V  V3  +  Pb  Vs 

4E  +  Pb  AV  +  Tb  Ss  (Tb,  Vs)  -  Tb  Sg  (Tb,  Vg)  =  S< 
„  1Hb 

s  “  Tv  +  Ss  (Tb'  V  -  Sg  (Tb'  V  ' 


With  the  energy  and  entropy  scales  correct t  the  dome  may  be 
constructed  by  solving  two  implicit  equations  for  selected  V 


P  (T,  V  )  =  P  (T,  V  ) 


Gs  IT,  Vs)  =  Gg  (T,  Vg)  , 


for  V  expressed  in  the  two  unknowns  T  and  Vr 

g  £ 


The  convergence  of  the  iteration  may  be  speeded  by  using  the  Clapyron 
equation  dP/dT  =  AH/TAV  to  estimate  the  temperature  T. 
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Calculation  of  Pressures 


Several  difficulties  can  be  foreseen  in  the  calculation  of 
pressure;  they  are  described  below. 

a.  Mixed-Phase  Region 

Clearly  there  are  two  possible  states  near  the  condensed 
vapor  line.  Consider  the  following  two  experiments  for  ambiguity. 

For  the  first  experiment,  put  the  solid  material  in  a  rigid 
container  with  a  piston  at  one  end,  not  attached  to  the  material. 
Heat  the  material  to  less  than  the  1  atm  melting  point,  and  then 
expand  the  volume  at  constant  temperature.  At  some  volume,  the 
piston  will  no  longer  touch  the  solid,  and  the  pressure  on  the 
piston  is  just  the  saturation  pressure  (or  vapor  pressure)  at 
that  temperature.  This  then  describes  a  material  under  the  liquid 
vapor  dome. 

The  second  experiment,  however  has  a  piston  that  is  rigidly 
attached  to  the  face  of  the  material.  As  the  piston  is  "expanded, 
the  material  will  go  into  tension  until  it  spalls.  Then  the  pres¬ 
sure  increases  rapidly  to  the  proper  point  in  the  liquid-vapor 
dome,  and  expansion  continues  as  in  the  first  experiment. 

With  these  two  experiments  in  mind,  then,  the  "mixed-phase" 

logic  in  a  hydro-code  becomes:  if  the  material  was  condensed, 

compute  pressure  as  solid  (P  ) .  If  P  is  less  than  the  spall 

s  s 

strength  (at  temperature  T) ,  compute  as  mixed  phase.  If  the 
material  was  a  gas,  or  if  one  surface  of  the  zone  is  free,  compute 
as  mixed  phase. 
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b.  Gas-Phase  Region  (V  >  Vcrit,  T  >  Tdome  (V)) 

In  this  region  there  is  no  difficulty,  provided  the  negative 
B  coefficient  (the  viral  coefficient)  does  not  dominate  the  pres¬ 
sure  equation  at  the  input  volume.  One  determines  that  the  sub¬ 
stance  is  a  gas  by  the  condition 

T-  a  (V.  )  ,  V. 
dome  in  m 

and  computes  the  pressure  by  solving  the  implicit  equation 
E  (T,  V.  )  =  E.  for  the  temperature.  The  pressure,  P,  is  given 
b?  P  =  p“<*.  4). 

c.  High-Density  Region  (V  <  Vcrifc,T  >  Tdome  (V)) 

One  problem  in  this  region  is  the  description  of  the  "transi¬ 
tion"  from  a  solid  to  a  high-density  gas.  Both  this  problem,  and 
one  other — that  there  is  no  guarantee  that  the  liquid-vapor  dome 
will  close  at  a  critical  point — can  be  traced  to  the  parameters  in 
the  solid  equation  of  state. 

If  the  liquid-vapor  dome  closes,  certainly  the  intersection 
of  the  solid  and  gas  lines  is  a  place  where  the  Helmholtz  free 
energies  are  equal,  and  continuity  arguments  would  say  that  the 
line  of  equal  Helmholtz  free  energies,  the  Tgub  line,  can  be  drawn 
from  vcr£t/  at  least  part  way  to  VQ.  Extrapolating  this  line, 
making  sure  that  E^  >  Eg  would  give  a  reasonable  transition  line. 
Remember,  however,  that  pressures  computed  near  this  artificial 
transition  line  may  not  be  accurate. 

If  the  dome  does  not  close,  there  is  no  assurance  that  the 
Tsub  can  ke  drawn  at  all.  There  are  two  ways  in  which  non¬ 

closure  may  occur.  The  first,  in  which  there  is  a  coexistence 
region  extending  to  high  temperatures,  is  shown  in  Figure  5. 
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0  V 

0  Volume 

FIGURE  5,  TWO -PHASE  DIAGRAM  FOR  SUBSTANCE  WITH 
NO  CRITICAL  TEMPERATURE  OR  VOLUME 


Figure  5  may  well  represent  the  condition  of  equilibrium 
between  a  solid  and  its  vapor — there  may  be  no  critical  point. 
This  is  certainly  reasonable  behavior,  and  may  be  pictured 
crudely  as  a  combination  of  the  vapor  pressure  increasing 
rapidly  with  temperature,  and  the  stiffening  of  a  material  at 
high  pressures  beginning  to  overcome  the  softening  associated 
with  high  temperatures. 

The  second  way  in  which  nonclosure  occurs,  however,  is 
the  case  where  there  is  no  temperature  at  which  the  Gibbs  free 
energies  equal,  shown  in  Figure  6. 


FIGURE  6,  TWO -PHASE  DIAGRAM  FOR  SUBSTANCE  WITH  CRITICAL 
TEMPERATURE  BUT  NO  CRITICAL  VOLUME 
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The  second  case  (shown  in  Figure  6)  is  more  difficult  to 
interpret  physically,  but  the  implication  is  that  since  the  two 
phases  may  not  coexist  above  a  certain  temperature,  there  must 
be  a  third  (separate)  phase  there.  Since  it  is  assumed  at  the 
outset  that  there  are  but  two  phases,  it  is  clear  that  the 
difficulty  is  with  the  separate  equations  of  state. 


In  the  first  case,  there  is  no  real  problem  in  calculating 

pressure.  There  is  no  transition  between  solid  and  vapor  phases 

in  the  high-density  region,  only  the  mixed  phase.  In  the  second 

case,  a  T^^  line  is  sought  starting  at  the  maximum  solid  volume 

by  the  standard  MONK  procedure.  If  there  is  no  solution,  a 

temperature  on  the  V  line  where  E  2.  E  is  used.  From  that  point 

S  Cf  s 

an  extrapolation  is  made  on  a  Tsub  line  at  constant  temperature 
(if  possible)  keeping  E^  >  Eg. 


d.  Alternative  Methods 

Certainly  one  alternative  to  the  methods  outlined  above  would 
be  tc  use  a  solid  equation-of-state  form  that  guarantees  the 
proper  behavior  at  the  critical  point.  This  was  the  method  taken 
in  Reference  15*  The  difficulty  with  (and  simplicity  of)  this 
method  is  that  there  are  no  adjustable  parameters  that  can  be  used 
to  accommodate  shock  wave  data.  It  would  seem,  then,  that  this 
approach  is  not  fruitful  for  use  in  a  hydro-code. 


e.  Consistency  Check 

After  the  boundary  of  the  dome  has  been  calculated,  one  severe 
test  of  consistency  is  that  the  critical  volume  and  temperature 
are  related  to  the  co-volume,  Bq,  and  temperature  e/k: 
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c/  k 

=  0.77 

V  ..t 

B 

o 

=  0.75 

V  •  *. 
cnt 

which  are  initially  estimated  to  scale  the  gas  equation  of  state. 
If  these  are  very  different  from  the  values  used  previously,  the 
dome  should  be  redrawn,  using  revised  values  of  the  co-volume  and 
characteristic  temperature. 


Major  areas  of  doubtful  values  of  pressure  using  this  formu¬ 
lation  would  be  near  the  transition  line,  T  ,  ,  and  near  the 

sub 

critical  point,  particularly  if  the  liquid-vapor  dome  does  not  close 

with  zero  slope.  It  should,  however,  do  well  with  condensation  and 

boiling  phenomena  for  P  <  and  for  release  adiabats  with 

P  >  P  ...  at  V  =  V 

cnt  cnt 


f .  Approximate  Version 

An  approximation  to  the  above  description  can  be  made  as 
follows.  Inside  the  liquid-vapor  dome,  the  two  relations  that 
follow  hold  rigorously. 


V 


+  E 


E 


(V  -  V  ) 
g  s 


V 


=  V 


(20) 


Figure  7  shows  the  reduced  vapor  pressure  as  a  function  of 
reduced  temperature.  This  is,  for  aluminum,  using  the  critical 
constants  T  =  6842°K,  P  =  4.7  kbar,  and  the  vapor  pressure  form 

C  w 

of  Reidel  (Reference  15)  > 

?v  --  exp  | a£n  t  t  C[36/t  -  35  -  t6  +  42  £n  t]} 
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In  further  analysis,  this  will  be  given  by  the  power  law 


P  =  t 
v 


6.85 


which  holds  reasonably  well  down  to  10  bar  (the  line  in  Figure  7) . 


Assuming  the  gas  to  be  ideal,  the  condition 


T  = 


E  -  E 
_g _ o 

3/2  R 


holds ,  where  E  is  the  internal  energy  of  the  gas  at  zero  tempera¬ 
ture  and  may  be  calculated  by  a  known  boiling  point  at,  for  example, 
1  atm. 


The  condition 


pv  («  - 


RT/V 


serves  to  determine  the  gas  volume  at  the  "dew  point."  If  we 
assume  further  that 

T  =  E  /3R  ( 

s 

i.e.,  small  volume  changes  in  the  solid,  T  ~  20,  we  may  make  the 
further  approximation  that 


V  -V 
g  s 


(v‘>[W¥  '"H 


n\V>  1  w/f  P/vii  /  O  A  \  — *  n  y>  <%V,4-  rt  ■» 

iJV^UU  UXUli  \4UV/  4.U4.  U41  l|g  /  wU  VMW.U4.li 


ng  E  -E 
g  s 
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Conditions  o£  Equations  (22)  and  (24)  give 


(E  -E  ) 

Es  =  3E  -3^-5-  =  2  <  W 


Defining  u  =  E  -  e  ,  E  =  2U,  Equation  (26)  becomes 


E-2U 


_______  -  E-2U 

"g  U+E-2U  eTu 
o  o 

Using  Equations  (21)  and  (23)  and  solving  for  V/V 

g 


(27) 


V. 


p  p 
c  v 

RT 


Pc  (t) 
RT 


6.85 


Substituting  Equation  (22)  for  T,  and  defining  U  =  3/2  RT 


V 

_  c  3/2  R 

V 

(  U 

V 

g 

R  U 

^3/2  RTc 

=  p  3V 

/5L 

,6.85 

c  2  U  - 

u_ 

) 

6.8! 


(28) 


Substituting  Equation  (28)  into  (25) 


Am6-85  3  v 

°\V  20 


and  equating  to 

_  E-2U 
Eq-U 


(27) 


which  on  rearranging  becomes 


P 


c 


/u_\6.85 

Iv 


2  U  |e-2u) 

3  V|E0^U) 


(29) 


Let  the  value  of  U  which  satisfies  (29)  be  U*.  Ti.en  in  the  mixed 
phase  region 
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T 


U* (E,V) 
"1/2  R 


1  11  /  n*  3 

\T  -  PC2U*^U^)  =  2U*  I,HS  EcIuat^on  (29) 


P  =  —  =  LHS  Equation  (29). 

g 


To  relate  this  to  a  hydro-code  calculation,  let  the  solid 
pressure  be  the  standard  polynomial,  P  ,  and  determine  the  approxi¬ 
mate  mixed-phase  region  by  Pg  (E,V)  <  Pv  (E/3  RTc)  and  if  V/Vg  <  1 
then  use  the  mixed-phase  logic. 


It  should  be  noted  that  the  approximations  restrict  the  range 

of  validity  to  T  <  T  .  ,  so  that,  for  aluminum  E.  £  1500  cal/g. 
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F.  CONCLUSIONS 


With  these  calculational  techniques,  one  can  predict  impulse 
m  materials  over  a  wide  range  of  deposited  energies.  For  high 
doses,  so  that  the  impulse  is  vapor  dominated,  tha  two-phase 
equations  of  state  in  hydro-codes  and  the  BBAY  expression  are 
reasonably  effective.  For  low  doses,  the  spall-dominated  region, 
the  BBAY  formula  is  inappropriate.  It  will  be  shown  in  Section  III 
that  the  Whitener  formula  with  an  appropriate  model  for  Ef  is 
effective  in  predicting  .impulse  when  proper  account  is  made  of  the 
deposition  time.  In  addition,  one-phase  equations  of  state  predict 
impulse  and  the  transmitted  stress  pulse. 

Equations  of  state  in  current  use  were  found  not  to  account 
for  the  effects  of  vapor-liquid  equilibria.  A  modification  and 
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extension  of  the  existing  McCloskey-Thompson  model  of  vapor-liquid 
equilibrium  effects  used  in  BBAY  integral  expressions  was  presented. 
Finally,-  a  technique  for  generating  a  thermodynamically  consistent 
model  for  vapor-liquid  equilibria  was  derived,  and  a  simplified 
version  suitable  for  direct  use  in  hydrodynamic  codes  was  shown. 

To  select  the  appropriate  models  for  predicting  impulse,  one 
needs  a  basis  of  empirical  information,  and  the  experiments  des¬ 
cribed  in  the  following  sections  were  performed  at  Physics  inter¬ 
national  on  the  738  Pulserad  to  provide  such  data.  Since  accurate 
electron  energy  deposition  profiles  are  required  for  all  of  tne 
models,  diagnosis  of  the  electron  energy  spectrum,  current,  and 
mean  angle  of  incidence  are  necessary  for  the  electron  transport 
(Monte  Carlo)  codes  used  to  calculate  the  deposition  profile. 

Section  III  presents  the  diagnostic  techniques  used  for  obtaining 
this  information,  as  well  as  the  techniques  for  measuring  impulse. 


41 


SECTION  III 


EXPERIMENTAL  TECHNIQUES 


A.  INTRODUCTION 


The  first  phase  of  the  experimental  portion  of  this  program 
consisted  of:  (1)  demonstrating  a  system  of  electron-beam  diag¬ 

nostics  to  determine  time-dependent  energy  deposition  in  the 
test  material  as  a  function  of  depth,  and  (2)  demonstrating  a 
reliable  technique  to  measure  electron-beam-generated  impulse. 

The  methods  used  to  accomplish  these  two  goals  are  discussed  in 
this  section.  The  beam  diagnostics  include  diode  voltage  and 
current  monitors  to  define  the  electron  beam  at  injection  into 
the  drift  chamber,  as  well  as  Rogowski  coils  and  a  Faraday  cup 
to  measure  the  net  and  primary  beam  currents  at  the  sample  loca¬ 
tion.  For  beam  intensities  sufficiently  low  to  preclude  material 
spallation  or  vaporization,  graphite  calorimeters  and  thin  foil- 
aluminum  dosimeters  were  used  to  measure  electron-beam  intensity 
and  in-depth  energy  deposition. 

The  ability  to  determine  energy  deposition  as  a  function  of 
material  depth  utilizing  voltage  and  current  monitor  outputs  in 
conjunction  with  mean  angles  of  incidence  of  beam  electrons  was 
independently  verified  by  the  following  comparisons: 

1.  Calculated  energy  deposition  with  that  measured  in  thin- 
foil  dosimeters. 

2.  Measured  electron  number  transmission  with  that  derived 
from  calculated  deposition  profiles. 

3.  Measure  rear-surface  pressure  in  a  material  with  a  well- 
known  equation  of  state  with  that  from  calculated  energy 
deposition  and  beam  intensity. 


Preceding  Page  Blank 
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lio  <mido  instrumented  with  Rogowski  coils  were 

employed  to  frail. .port  electron  beams.  The  geometry  of  the 

guide  was  designed  to  shape  the  beam  intensity  so  that  the 

tlueno-  ical/cm1')  is  uniform  over  the  guide-cone  exit  area 

2 

(approximately  3.U  cnt“i  .  The  Rogowski-coil  instrumentation 
provides  a  measurement  of  the  net  current  of  the  beam,  which 
can  be  empirically  related  to  the  primary  current.  This  re¬ 
lationship  has  been  used  to  predict  beam  intensity  (fluenco) 
within  ±10%. 

A  ballistic  pendulum  was  used  to  measure  electron-beam- 
generated  impulse  in  test  materials.  The  important  features 
of  the  pendulum  were  a  variable  transformer  recording  the 
time  history  of  the  pendulum  deflection  and  a  low  damping  co¬ 
efficient  of  oscillation.  The  guide  cone  and  pendulum  system 
are  mounted  on  an  inertial  platform  used  to  decouple  the 
measuring  apparatus  from  "bulkhead  shock"  caused  by  the  shock 
generated  in  the  oil  switch  between  the  transmission  line  and 
the  diode.  The  beam  guides  have  been  slotted  to  minimize  the 
effects  of  anode  debris  on  pendulum  motion. 

B.  DIAGNOSTIC  TECHNIQUES 

Electron-beam  parameters  are  monitored  in  two  regions: 
the  diode  and  the  drift  chamber.  Diode  diagnostics  provide 
a  time  history  of  electron  accelerating  voltage  and  current 
at  the  anode  plane.  This  information  is  required  to  determine 
the  characteristics  of  beam  production  in  the  diode. 

Following  injection  of  the  beam  into  the  drift  chamber, 
the  characteristics  of  beam  transport  to  the  sample  location 
can  alter  the  energy  spectrum  generated  in  the  diode.  Addi- 
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tional  diagnostics  in  the  drift  chamber  are,  therefore,  required 
to  specify  the  electron-beam  environment  at  the  sample  location. 
Strong  evidence  indicates  that  the  voltage  waveform  is  noL  signi¬ 
ficantly  altered  for  beam  transport  lengths  up  to  30  cm.  However, 
the  primary  current-pulse  shape  can  undergo  substantial  change 
caused  -by  losses  of  electrons  with  high  transverse  energy  com¬ 
ponents  as  well  as  by  the  "erosion"  of  the  beam  front  (Reference  3) . 

What  follows  is  a  description  of  the  monitors  in  the  diode 
and  drift  chamber  regions,  including  their  calibration  and  function 
to  determine  electron-beam  characteristics . 

1.  Diode  Diagnostics 

The  Physics  International  738  Pulserad  was  used  to  generate 
the  high  energy  density  states  required  in  this  program.  The 
pulser  consists  of  a  38-stage  Marx  generator,  used  to  pulse-charge 
an  8.5-ohm,  oil-filled  coaxial  transmission  line.  The  trans¬ 
mission  line  is  switched  into  the  load,  a  field  emission  diode 
consisting  of  a  2^-in.  diam,  600-needle  cathode,  and  a  thin, 

0 .00025-in. -thick,  aluminized  Mylar  anode.  Diode  impedance  is 
variable  between  «  1  and  «  8  ohms  by  adjustment  of  the  anode-to- 
cathode  spacing,  such  that  at  typical  pulse  charge  voltages  of 
3.5  mV,  the  machine  output  can  be  varied  in  the  2C0  keV  to  1  MeV 
mean  electron  energy  range  at  current  levels  of  250  to  125  kA, 
respectively. 

The  diode  diagnostics  consist  of  a  resistive  voltage  divider 
attached  to  the  cathode,  a  self-integrating  loop  current  monitor, 
and  an  open  loop  dl/dt  monitor.  These  three  monitors  are  situated 
in  the  diode  (Figure  8) .  Their  construction  is  described  in 
Reference  2. 
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FIGURE  8,  ANODE-CATHODE  REGION  AND  PLACEMENT  OF  DIAGNOSTICS 
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These  monitors  were  calibrated  as  follows: 


1.  The  diode  current  monitor  was  calibrated  by  firing  the 
beam  into  a  Faraday  cup  located  at  the  anodo  plane  at 
both  high  (200  kA)  and  low  (50kA)  current  levels.  The 
diode  voltage  monitor  was  calibrated  at  25  kV  with  an 
external  voltage  source. 

2.  The  inductive  pickup  of  the  voltage  monitor  (Figure  9c) 
was  measured  for  anode-cathode  shorted  shots  and  cor¬ 
related  with  the  signal  from  the  dl/dt  probe  located  in 
the  diode  region.  The  dl/dt  probe  was  then  used  as  a 
measure  of  the  inductive  component  of  the  diode  voltage 
signal  on  all  subsequent  shots.  This  component  was  sub¬ 
tracted  from  the  voltage  monitor  signal  (Figure  9 d )  to 
yield  the  actual  accelerating  voltage  pulse  shape 
(Figure  9e) . 

3.  A  series  of  shots  was  fired  into  total  shopping  graphite 

calorimeters  near  the  anode.  The  total  beam  calories 

deposited  in  the  calorimeter  agreed  (±5%)  with  /”  VI  dt 

0 

calculated  from  diode  voltage  and  current  records . 

2 .  Beam  Transport 

The  beam  in  the  drift  chamber  is  controlled  by  metallic  guide 
cones  (Figure  10) .  These  cones  produce  an  extremely  flat  fluence 
distribution  over  the  3/4-in.  exit  diameter  (at  drift  chamber 
pressures  in  the  range  0.5  to  0.75  Torr)  as  confirmed  by  uniform 
depth  cratering  in  aluminum,  and  X-ray  pinhole  photography  of  the 
bremsstrahlung  produced  in  a  tantalum  target  at  the  cone  exit. 

The  pinhole  photographs  in  Figure  11  show  uniform  density  spots 
corresponding  to  the  area  at  the  cone  exit.  The  elliptical  shape 
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FIGURE  9.  DIODE  DIAGNOSTIC  TRACES 

(TIME  SCALE  =  20  nsec/cm) 
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FIGURE  10.  SLOTTED  ELECTRON  BEAM  GUIDE  CONE  INSTRUMENTED 
WITH  TWO  ROGOWSKI  COILS 


of  the  spots  results  from  positioning  the  pinhole  camera  30 
degrees  off-axis  from  r.he  center  line  of  the  cone.  Figure  12 
shows  the  top  and  side  view  of  an  electron-beam-caused  crater 
in  an  aluminum  Blab.  The  uniform  crater  depth  confirms  a 
uniform  fluence  distribution  over  the  exit  area  of  the  cone. 


The  level  of  the  uniform  fluence  distribution  (at  fluence 
levels  low  enough  tv.  preclude  removal  of  graphite  by  the  beam) 
is  determined  by  allowing  the  beam  to  impinge  on  a  totally 
absorbing  (5/16-in.  thick)  segmented  calorimeter  array  placed 
at  the  cone  exit.  The  energy  stopped  in  those  blocks  which 
are  well  within  the  projection  of  the  exit  area  of  the  cone 
onto  the  calorimeter  array  (Figure  13)  is  divided  by  the  frontal 
area  of  the  blocks  ( H  x  %-cm  )  to  yield  the  fluence  level  (cal/cm  ) 

The  fluence  level  indicated  by  blocks  near  the  circumference 
of  the  cone  exit  area  projection  is  underestimated  for  two  rea¬ 
sons.  These  blocks  are  not  irradiated  over  their  entire  area 
(fluence  calculations  average  incident  calories  over  entire  front 
surface  of  the  block)  and  also  they  lie  on  the  slope  of  a  severe 
temperature  gradxent  from  the  hotter  blocks  in  the  center  to  the 
colder  blocks  near  the  edge  of  the  calorimeter  array.  The  heat 
transfer  is  demonstrated  by  fluences  recorded  in  blocks  that  were 
not  irradiated  (see  Figure  13,  Calorimeter  Array) .  Consequently, 
the  fluence  map  obtained  from  the  calorimeter  array  is  a  poor 
indication  of  fluence  uniformity. 


Figure 
calorimetry 
portions  of 


13  shows  a  radial  fluence  distribution  obtained  from 
incorporating  a  correction  for  irradiation  cf  only 
blocks  along  the  cone  circumference  projection. 
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Fluence  (cal/cm 


3. 


D riffcCh antbe r  Di-ag n o sties 


The  beam  guide  cones  are  instrumented  with  two  Rogowski 
coils,  each  of  which  is  a  loosely  wound  torus.  The  Rogowski 
coil  is  self-integraiing  and  has  a  voltage  output  proportional 
to  the  average  flux  threading  all  its  turns.  This  diagnostic 
measures  tne  net  current  flowing  through  the  guide  cone.  A 
completed  coil  is  shown  in  Figure  14.  Neglecting  any  mutual 
inductance  between  the  separate  turns  of  the  Rogowski  coil, 
the  output  voltage  is 

V  =  ~  nAp  *s7ra  (mks) 

L  o 

where  R  =  resistance  of  integrating  resistor, 

L  =  net  self-inductance  of  the  coil, 
n  -  number  of  turns  in  coil, 

A  =  area  of  a  single  turn, 
a  =  mean  radius  of  ehe  torus. 

In  view  of  the  approximations  involved  in  determining  the 
Rogowski  coil  output,  an  empirical  calibration  was  necessary. 

A  calibration  fixture  and  fast  readout  system  have  been  developed 
from  the  following : 

1.  Fast  risetime  pulse  generator,  10  A  into  50-ohm  load 
with  a  2-nsec  risetime. 

2.  Sampling  oscilloscope,  2  mV/cm  sensitivity,  0.3-nsec 
risetime. 

3.  Impedance  matching  load  threading  the  coil. 

A  typical  calibration  trace  on  the  sampling  oscilloscope  is 
shown  in  Figure  15a,  where  the  upper  trace  is  the  coil  output,  and 
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the  lower  is  the  input  current  pulse.  Figure  15b  shows  the  coil 
response  to  a  longer  (270  nsec)  10-A  input  pulse.  A  summary  of 
the  characteristics  for  coil  15  is  given  below. 


COIL  CHARACTERISTICS 

2  cm  i.d.  Rogowski  Coil  Number  15 
Risetime  0.9  nsec 
Absolute  Sensitivity  4.8  V/kA 
Decay  Constant  409  nsec 


The  aim  of  the  inclusion  of  Rogowski  coils  in  the  beam  guide 
cones  was  the  prediction  of  fluence  levels  when  a  sample  completely 
covers  the  cone  exit  area  thus  making  fluence  data  collection  im¬ 
possible.  A  Faraday  cup  is  used  to  achieve  this  predictive  capa¬ 
bility  t Figure  16) .  The  Faraday  cup  is  a  charge  collector  combined 
with  an  integral  millj.ohm  shunting  resistor  and  incorporates  a  fast 
vacuum  pumping  system  and  a  conical  collector  specifically  designed 
to  intercept  high-current,  low-energy  electron  beams.  A  thin  rub¬ 
ber  membrane  acts  as  a  vacuum  seal  in  front  of  the  collector,  and 
in  principle,  prevents  the  plasma  formed  in  the  beam  channel  from 
shorting  out  the  current  path  through  the  monitoring  resistor. 
Typical  output  signals  of  the  Rogowski  coil  and  the  Faraday  cup 
are  shown  in  Figure  17. 

With  1^,  1 2,  and  I™ax  defined  in  Figure  17,  the  following 
empirical  relation  between  the  primary  and  net  current  for  a  fixed 
drift  chamber  pressure  {e.g.,  K  =  2.0  at  0.5  Torr)  has  been  ob¬ 
served  : 


I 
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Rcgowski  Coil  Output  Parameters 


The  rms  deviation  in  the  prediction  of  the  primary  current  peak 
amplitude  is  tf> 1  (Figure  13). 

Since  a  Faraday  cup  trace  of  the  primary  current  is  not 
available  during  a  sample  irradiation  shot,  the  time  dependence 
of  the  primary  current  was  approximated  by 

tp{t)  =  Ic  sin  (t  -  a1) 

where  1  =  (K82  ImaX 

o  p 

c'  =  full  width  -  86  nsec 

a'  =  lag  time  compared  with  V(t  =  0)  =5  nsec 

The  accelerating  voltage,  V(t) ,  from  the  voltage  monitor  trace 
was  approximated  by  a  trapezoid: 

V, 

V  ( t )  =  t  for  0  s  t  g  a  =  12  nsec 

a 

(V  -  V,)  t  -  (aV  -  bV  ) 

V(t'  «  — - t- - for  a  <;  t  £  b  =  50  nsec 

a  -  a 

V  (c  -  t) 

V  ( t)  =  — (cPbl ^or  ^  £  t  £  c  =  86  nsec 

where  and  V2  are  defined  in  Figure  19: 


60 


A  k  -  1.59  at  0.15  Torr 


FIGURE  18  MEASURED  PEAK  PRIMARY  CURREN1 ,,  I  ma:<  VERSUS 
THE  I,  +  kl-  CALCULATION  P 


Using  the  approxima  :ions  for  V(t)  and  I  (t)  as  defined 

It 

abov';,  let 

CM 

3C(cal)  “  4.186  f0  V  (t)  Ip  (t)  dt  =  IQ(kA)  [3.12  V1  (kV)  +7.34  (kV)]x  10 

represent  the  total  beam  calories  near  the  cone  exit.  The  uniform 
fluence  near  the  cone  exit  now  permits  division  by  the  cone  exit 
cross-sectional  area  giving  an  approximate  fluence  level 

?(cal/cm2)  5  =  I  (kA)  [1.09  (kV)  +2.55  V2  (kV)  ]  x  10-3 

where  IQ  =  0.82  I™ax  =  0.82  (I1  +  2.0  I2) ,  with 

k  =  2 . 0  at  0 . 5  Torr 

A  =  2.85  cm2 

<jf  depends  only  or.  parameters  from  the  net  current  (Rogowski 
coil  output)  and  the  accelerating  voltage  (diode  voltage  monitor 
output) .  This  information  is  available  on  any  shot,  including 
sample  irradiation. 

The  results  of  a  series  of  fluence  shots  into  segmented 
graphite  calorimeters  were  used  to  test  the  accuracy  of  the  total 
beam  calories  prediction,  3C.  The  rms  deviation  from  the  measured 
total  beam  calories,  3C,  in  the  range  of  60  to  170  cal  was  ±10% 
(Figure  20).  A  further  comparison  consisted  of  measuring  the 
primary  beam  current  with  a  Faraday  cup  at  the  guide  cone  exit. 

The  waveforms  of  the  voltage  and  primary  current  were  then  used 
in  an  exact  integration  of  /V(t) I„(t)dt  ro  obtain  the  total  beam 

tr 

calories  which  were  compared  with  3C  calculated  from  approximations 
to  those  waveforms.  These  data  are  also  presented  in  Figure  20. 
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The  mean  electron  energy  is  defined  as 

CO 

I  V  (t.)  I  (t)dt 

*o 

-  “to 

/  I  (t )dt 
J o 


and  3C 


/ 


V  (t)  I  (t)dt  with  I(t)  =  Io  sin  -n/c  '  (t 


-  a') 


Taking  advantage  of  the  ^-approximation ,  the  mean  electron 
energy  can  be  estimated  by 


V (t) I (t)dt 


c 


00 


I (t)dt 


0.266  (kV)  +  0.626  V2  (kV)  <£> 


The  correlation  between  <E>  and  (e)  can  be  used  as  a  further 
check  on  the  validity  of  the  ^-approximation.  Figure  21  is  a  plot 
of  <E>  versus  (s),  when  <E>  was  calculated  from  the  signals  of  the 
voltage  and  current  monitors  without  approximation.  The  rms  de¬ 
viation  of  this  comparison  is  less  than  ±10%. 


C.  ELECTRON  ENERGY  DEPOSITION  PROFILE  AND  SUMMARY  OF  BEAM 
DIAGNOSTICS 

1 .  Depth  Dose  and  Electron  Number  Transmission  Measurements 

The  method  used  to  predict  the  deposition  profile  is  diagramed 
in  Figure  22.  The  electron  beam  is  fired  into  an  array  of  0.003-in. 
aluminum  foils  with  thermocouples  attached  (uepth  dose), to  record 
the  temperature  rise.  The  data  are  translated  into  energy  deposited 
in  each  foil,  and  a  normalization  is  obtained  by  dividing  the  cal/g 
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U)  kev 

FIGURE  21  .  EXACT  CALCULATION  OF  THE  MEAN  ELECTRON 
ENERGY,  <E> ,  VERSUS  THE  APPROXIMATION, 
(p\  =  0.266  V,  +  0.626  V. 


FIGURE  2 2.  DIAGNOSTICS  LOGIC  DIAGRAM 


deposited  in  each  foil  by  the  incident  fluence.  The  fluence  is 
defined  by  the  total  calories  stopped  in  the  entire  array  divided 
by  the  collimator  are  a.  The  result  is  a  measured  depoe.  ..tion  pro¬ 
file  for  aluminum  in  histogram  form. 

For  the  same  shot,  the  time-dependent  accelerating  voltage, 

V(t),  and  the  diode  monitor,  I(t)  (Reference  3),  are  divided  into 
small,  equal,  time  increments,  At,  generating  a  set  [V(t^),  I(t^)], 
which  becomes  the  input  to  an  electron  deposition  code  based  on  an 
interpolation  of  Spence's  data  (Reference  2).  The  code  considers 
each  pair,  [V(t^),  I(t^)],  to  be  a  monoenergetic  source  of  energy 
eV(t^)  and  the  electron  number  At  I(t.)  from  which  it  calculates 
the  spectrum  and  the  net  deposition  profile  in  aluminum.  The  de¬ 
position  profile  obtained  carries  the  normalization  of  unit  fluence, 
i.e.,  eal/g/eal/eni  .  The  depth-dose  histogram  is  then  compared  with 
the  deposition  profile  calculated  from  the  accelerating  voltage  and 
diode  monitor  traces.  Figures  23  and  24  show  the  calculated  de¬ 
position  profiles  for  mean  electron  energies  of  283  and  216  keV, 
respectively.  Superimposed  over  these  curves  are  the  measured 
depth  dose  histograms. 

A  consistent  trend  appears  in  the  comparison  of  the  deposition 
profiles.  The  first  foil  always  reads  a  value  higher  than,  and  the 
second  foil  a  value  lower  than  that  predicted  from  the  diode  diagnos¬ 
tics.  This  behavior  is  consistent  with  non-paraxial  electron  tra¬ 
jectories  and  appears  to  be  characteristic  of  all  high-current 
beams.  Clearly,  a  precise  definition  of  electron  trajectories  is 
required  to  resolve  this  question,  and  experiments  !;o  this  end 
have  been  conducted  under  DASA  Contract  No.  DASA-01-68-C-Q096 , 
Specifically,  these  measurements  incorporated  a  Faraday  cup 
behind  filters  of  various  thicknesses,  and  the  transmission  data 
were  compared  with  Monte  Carlo  electron  deposition  calculations 
(Figures  25  and  26). 


67 


al/ 


283  kcV 


50  nsec  100  nsec 


50  nsec  100  nsec 


Note: 


0  Monte  Carlo,  (9)  =0° 
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FIGURE 


ELECTRON  NUMBER  TRANSMISSION  THROUGH 
ALUMINUM  FILTERS 


0 


Kv.'S'i Its  of  these  investigations  indicate  a  good  correlation 
between  the  measured  mean  incidence  angles,  (0),  and  these  calcu¬ 
lated  from  the  balance  between  electron  transverse  and  magnetic 


pressures  (Reference  3) ,  which  predicts 


+  1-1 


_ Peak  Nc-t  Current 

Peak  Primary  Current 


Data  from  a  recent  program  funded  by  Sandia,  Livermore  (Reference  9), 
illustrate  the  correlation  between  measured  and  calculated  (0)  at 
the  sample  location: 

(E>  =  480  keV 

a-im)  =  o.6 

:net  ”  15  kA 
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the  data 


<bt2> 

<e> 


25  kA 


0.38 

31° 


The  calculated  angle  is  in  good  agreement  with 
shown  in  Figure  25. 


2.  Summary  of  Beam  Diagnostics 

The  techniques  of  determining  the  electron  fluence  level 
and  normalized  deposition  profile  as  discussed  in  the  preceding 
section  are  summarized  in  Figure  27.  A  brief  review  follows: 

The  scaled  dl/dt  probe  signal  representing  the  inductive 
pickup  of  the  voltage  monitor  is  subtracted  from  the  voltage 
monitor  signal,  yielding  the  accelerating  voltage  across  the 
anode-cathode  space.  This,  in  addition  to  diode  current  and 
the  mean  angle  of  incidence,  is  used  as  the  input  to  Monte  Carlo 
electron-deposition  calculations  yielding  the  mean  electron 
energy  and  a  deposition  profile  normalized  to  unit  flucnce. 

Using  the  signals  from  the  voltage  monitor  and  the  Rogowski  coil 
(net  beam  current  near  the  sample) ,  the  fluence  during  an  actual 
sample  irradiation  shot  is  determined  and  used  to  scale  the  nor¬ 
malized  deposition  profile  to  yield  the  actual  deposition  profile 
for  that  shot. 


Data  collected  in  a  recent  program  {Reference  5)  has  afforded 
an  excellent  opportunity  to  test  the  validity  of  the  738  Pulserad 
diagnostics  in  defining  the  electron-beam  environment.  Figure  28 
shows  the  comparison  between  calculated  and  measured  stress  his- 
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APPLICATION  OF  ELECTRON  BEAM  DIAGNOSTICS  TO  MATERIAL 


tori.es  and  Figure  29  gives  a  plot  of  peak  stress  (measured  and 
calculated)  versus  fluence  in  beryllium.  The  agreement  between 
the  measurements  and  calculations  indicates  a  good  confidence 
level  in  the  diagnostics. 

D.  IMPULSE  MEASUREMENT  TECHNIQUES  AND  DATA  SUMMARY 

The  ballistic  pendulum  used  in  this  program  is  shown  in 
Figure  30.  The  pendulum  axle  supporting  the  rod  and  bob  rotates 
on  two  instrumentation  bearings  and  is  mechanically  coupled  to 
a  variable  transformer  to  record  the  pendulum  deflection.  Ro¬ 
tation  of  the  pendulum  axle  changes  the  magnetic  coupling  be¬ 
tween  an  externally  powered  primary  and  a  secondary  coil  whose 
output  signal  is  then  displayed  on  an  oscilloscope.  Figure  31 
shows  a  typical  output  signal  recording  the  pendulum  deflection 
as  a  function  of  time.  The  variable  transformer  response  was 
calibrated  by  photographing  its  output  at  successive  pendulum 
deflections  in  five-degree  increments  (Figure  32). 


Figure  33  shows  the  calibration  curve  obtained  by  plotting 
pendulum  deflection,  9,  versus  he  output  signal  from  the 
variable  transformer.  The  output  is  linear  for  0°  s  8  s  20°, 
with  increasing  nor?  -linearity  for  0  >  20°.  Although  the 
calibration  curve  permits  meaningful  operation  in  the  U  >  20° 
region  of  maximum  deflection,  the  normal  operating  procedure  is 
to  increase  the  bob  mass  of  the  pendulum  in  advance  of  antici¬ 
pated  deflections  of  0  >  20°.  Low-impulse  sensitivity  is 
obtained  by  decreasing  the  bob  mass  while  increasing  the 
oscilloscope  sensitivity. 

The  differential  equation  governing  the  motion  of  a  damped 
simple  pendulum  is  of  the  form 

0  +  2B9  +  ^  sin  0=0 
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O  Quartz  Gauge  Data 
A  POD  Calculations 


l - — - 1 - 1 - - 1 — 

0  50  100  150 

2 

Fluence  (cal/cm  ) 

FIGURE  29.  PEAK  STRESS  VERSUS  FLUENCE  QUARTZ  DATA  FOR 
BERYLLIUM 


FIGURE  30  t  BALLISTIC  PENDULUM 


Upper : 

6.67  deg/cm 
100  msec/cm 

Lower : 

1.67  deg/ cm 

10  msec/cm 


Pendulum  Deflection.  Second  Peak  Is  A  Rebound  From  A  Bob  Stop 

To  Prevent  Sample  Impact  On  Guide  Cone 


FIGURE  31.  RECORDED  PENDULUM  DEFLECTION 


100  msec/cm 


FIGURE  32.  VARIABLE  TRANSFORMER  CALIBRATION  AT  A0 


FIGURE  33.  VARIABLE  TRANSFORMER  CALIBRATION  CURVE 
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where  the  damping  is  proportional  to  the  angular  velocity,  o. 

For  angles  up  to  20  degrees  (sin  0) / 0  «  0.98.  Thus,  to  a  good 
approximation 

0  +  260  +  |  0  =  0 

It  can  be  shown  that  0 ( t )  =  Ae  ^  sin  ut  with  i»  =  y/g/  z  -  p 
satisfies  the  equation  of  motion. 

A  determination  of  the  damping  coefficient,  6,  is  made  by 
recording  a  number  of  free  oscillations  of  the  pendulum  resulting 
from  electron-beam-generated  impulse.  Comparing  peak  deflections, 
a  number  of  oscillations,  n,  apart  yields 

-6t. 

6(t^}  Ae  sin  wt^ 

0 (t  )  =  -8tl 

Ae  sin  u)ty 

Thus,  the  damping  coefficient 

w  1  0(fcl) 

0  2nn  ln  0(t2)  nT  ln  0 (t2) 

where  the  period  T  and  the  ratio  0(t^)/O(t2)  are  taken  directly 
from  the  output  of  the  variable  transformer.  Damping  coefficient 
determinations  are  performed  prior  to  every  experimental  run.  A 
typical  value  is  8  =  0.03  sec  .  The  period,  T,  of  the  pendulum 
is  0.95  sec  for  most  bob  masses  used.  The  first  peak  deflection 
from  the  beam-generated  impulse  occurs  at  t  =  1/4T  =  0.24  sec. 

The  damping  factor,  under  these  conditions,  has  the  value 
e~p,t  =  0.99  «  1.0. 


-6 (t,  -  t„) 


=  e 


e2irn6/w 


82 


Consequently,  bearing  drag  associated  with  pendulum  motion, 
although  small,  can  be  accounted  for  in  the  impulse  calculation. 


Applying  the  initial  condition  SdO/dt  =  v  at  t  =  0  to  the 
solution  of  the  damped  pendulum  equation  of  motion  yields 


G(t)  =  —  {«  +  B^/w)  e  ^  sin  wt 
9 


But  since  typically 


„2  _  (0.03;  (0.95) 

6  /<o - - - 


=  0.0014  <<  to, 


0(t)  =  —  e  sin  wt 
9 


2vv  -(St 

gijr-  e  sm  wt 


The  maximum  deflection,  9  ,  occurs  at  t  =  1/4  T  where  sin  wt  =  1. 

IiImX 

Solving  the  0  equation  for  v  and  multiplying  by  the  total  moving 
max 

mass  of  the  pendulum,  the  momentum  imparced  to  the  pendulum  is 


Mg  T  0  e  8(T/4) 

(MV) - !!§2 -  =  2.72  MT  e^Mc-g)  exp  <|2) 

The  peak  deflection,  0  ,  and  the  period  of  the  pendulum,  T,  are 

max 

obtained  directly  from  the  oscilloscope  trace  of  the  output  from 
the  variable  transformer  that  monitors  the  pendulum  motion 
(Figure  31) . 


Consideration  of  the  initial  slope  of  the  pendulum  deflection 
(Figure  31,  lower  trace)  affords  another  calculation  of  the  impulse 
impacted  to  the  pendulum,  if  i  is  the  pendulum  center  of  mass, 
then 


A (MV)  =  M& 


A0 

At 


2 


w 


+  8 


AG 
2  It 


Mg  A0 
2  At 


w 
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Since  the  damping  coefficient  measurements  cannot  reflect 
effects  from  static  frictional  forces  at  the  start  of  the 
pendulum  swing,  a  calibration  of  the  pendulum  with  a  known 
impulse  was  necessary.  A  CC^-cartridge-fed  air  pistol  provided 
the  calibrating  impulse.  The  pellet  velocity  was  obtained  from 
the  transit  time  between  two  photocells  separated  by  a  known 
distance;  the  pellet,  mass  was  measured  with  an  analytical  balance. 
The  pellet  was  brought  to  rest  in  the  pendulum  bob  by  a  layer  of 
absorbing  material.  The  setup  is  shown  in  Figure  34  and  the 
calibration  data  are  shown  in  Figure  35. 

Because  of  the  motion  of  the  738  Pulserad  during  an  electron 
beam  pulse,  it  became  necessary  to  mount  the  ballistic  pendulum 
and  beam  guide  cone  on  a  stable  platfon.  The  experimental 
setup  is  shown  in  Figure  36. 

Ex1'  sting  pendulum  bobs  accommodate  disk-shaped  samples  of 
1.30  in  diam  and  thicknesses  from  0.080  to  0.25  in.  This, 
however,  is  not  a  limitation,  since  bobs  are  easily  fabricated 
to  suit  the  task.  The  samples  are  clamped  between  two  retaining 
rings  and  this  essentially  constitutes  the  front  face  of  the 
pendulum  bob  (Figure  37) . 


The  accuracy  of  impulse  measurements  depends  on 
quantities:  the  momentum  imparted  to  the  bob,  A (MV) 
crater  area  of  the  sample,  A.  The  error  in  A (MV)  is 
demonstrated  in  the  calibration  (Figure  35)  iri  which 


two 

and  the 
best 

the  deviation 


34 


Variable 

Transformer 


Air  Blast  Shield 


Enoto  Cells 
With  Light  Sources 


-  Light 

..  l/Q-’in.  Aperture 
•Photocell 


+  12 


Systron-Donner  Frequency  Counter 
(Model  1034) 

Trigger  for  Oscilloscope  Monitoring 
Variable  Transformer 


FIGURE  34  PENDULUM  CALIBRATION  SETUP 


FIGURE  35  GUN  PELLET  MOMENTUM,  mv  VERSUS  MOMENTUM 
RECORDED  BY  PENDULUM  L  (MV) 


Diode  Current 
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Variable 
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FIGURE  36.  IMPULSE  MEASUREMENT  SETUP 


from  the  calibrating  momentum  is  less  that  ±  3%.  This  error  bar 
includes  inaccuracies  in  the  measurement  of  the  total  moving  mass 
and  the  error  incurred  in  measuring  the  maximum  amplitude  and 
period  of  oscillation  from  the  oscilloscope  photograph  of  the 
variable  transformer  output.  The  error  in  the  measurement  of 
beam  crater  area  can  be  minimized  by  photographing  samples  "head- 
on"  and  measuring  A  from  the  pictures  using  a  polar  planimeter. 

The  deflection  readout  of  the  pendulum  is  unaffected  by 
electron-beam-generated  noise  since  the  bob  remains  stationary 
during  deposition  time  and  begins  to  move  long  after  all  noise 
has  subsided.  Variations  in  the  deposition  profile  shape  affect 
the.  pendulum  only  through  the  magnitude  of  the  total  impulse  they 
generate.  The  pendulum  deflection  can  be  controlled  by  varying 
the  mass  of  the  bob.  Oscillations  on  the  order  of  a  tenth  of  a 
degree  are  easily  measured  by  increasing  oscilloscope  sensitivity 
to  the  variable  transformer  signal.  There  is  virtually  no  upper 
limit  to  measurable  impulse  since  bob  mass  can  be  increased  to 
the  kilogram  range. 

A  pendulum  with  a  manganin  gauge  built  into  the  bob  is 
currently  under  development.  Damping  characteristics  like  those 
of  the  standard  pendulum  have  been  achieved.  The  remaining 
difficulty  is  shielding  the  gauge  from  the  RF  environment  in  the 
electron  beam  drift  chamber  without  coupling  the  machine  motion 
to  the  pendulum.  This  problem  is  presently  being  addressed  and 
once  solved  will  allow  the  measurement  of  total  impulse  and  rear 
surface  stress  simultaneously. 
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The  techniques  described  above  have  been  used  to  generate 
impulse  data  in  aluminum.  These  data  are  summarized  in 
Table  II.  Since  the  data  cover  a  range  of  mean  energies,  a 
number  of  typical  energy  deposition  profiles,  normalized  to 
unit  fluence,  are  shown  in  Figures  38  to  43. 

The  maximum  front  surface  dose  in  these  experiments  was 
1800  cal/g.  This  is  shown  in  the  following  section  to  be 
within  the  range  of  the  spall-dominated  impulse.  In  Section 
IV,  a  model  is  presented  for  this  regime,  which  is  in  good 
agreement  with  the  data.  In  addition,  measurements  and  cal¬ 
culated  values  of  transmitted  stress  and  ejected  velocities 
are  presented. 
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TABLE  II 


ALUMINUM  IMPULSE  DATA 


<E> 

(keV) 

<$> 

(cal/cm2) 

A2 
(cin  ) 

A  (MV) 

(KILO-DYNE-SEC) 

I 

(ktap) 

Am 

(g) 

(208) 1 

292 

2.73 

4.08 

1.51 

0.082 

(222) 

(30) 

3.2 

5,93 

1.85 

0.109 

220 

31 

3.0 

5.45 

1.83 

0.119 

(225) 

35 

3.2 

5.55 

1.72 

0.131 

(232) 

(37) 

3.3 

7.66 

2.32 

0.151 

(233) 

41 

3.35 

7.19 

2.15 

0.117 

(260) 

41 

3.0 

6.54 

2.18 

0.140 

219 

42 

3.5 

7.19 

2.07 

0.138 

249 

43 

3.2 

6.02 

1.88 

0.114 

257 

43 

2.8 

7.30 

2.61 

0.154 

(223) 

45 

3.3 

6.76 

2.05 

0.137 

353 

52 

(~4 . 0 ) 

13.17 

~  3 . 3 

0.242 

266 

60 

3.3 

9.05 

2.74 

0.155 

371 

60 

3.9 

17.08 

4.40 

0.279 

(386) 

61 

(~3.8) 

13.63 

~  3 . 6 

0.241 

620 

113 

3.7 

31.25 

8.44 

0.537 

Notes  to  Tabular  Data 


Column  1, 

Column  2. 
Column  3. 


Parentheses  indicate  mean  energy  estimated  by  approxi¬ 
mate  integration  scheme,  others  by  numerical  integration 
of  voltage  and  current  traces. 

Parentheses  indicate  fluence  estimated  by  approximate 
integration  of  current  trace. 

Parentheses  indicate  uncertain  measurements  of  area  due 
to  irregular  crater  edge. 
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Deposited  Energy 
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FIGURE  38.  NORMALIZED  DEPOSITION  PROFILE 
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NORMALIZED  DEPOSITION  PROFILE 
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FIGURE  40.  NORMALIZED  DEPOSITION  PROFILE 
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Mean  Angle  =  52.9°  from  Normal 
Mean  Electron  Energy  =  0.434  MeV 


Distance  (cm) 


FIGURE  41.  NORMALIZED  DEPOSITION  PROFILE 
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FIGURE  42.  NORMALIZED  DEPOSITION  PROFILE 


Mean  Angle  =  58.6  from  Normal 
Mean  Electron  Energy  =  0.620  MeV 
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FIGURE  43.  NORMALIZED  DEPOSITION  PROFILE 


SECTION  IV 


EXPERIMENTS  AND  ANALYSIS  ON 
MELT-DOMINATED  IMPULSE 


A.  INTRODUCTION 

The  two  previous  sections  describe  analytical  and  experi¬ 
mental  techniques  used  for  studying  impulse  generation  and  the 
impulse  data  that  were  collected  for  aluminum  in  the  melt- 
dominated  or  liquid  splash-off  regime. 


In  this  section,  a  model  for  melt-dominated  impulse  is 
developed  that  gives  good  agreement  with  the  data  collected  here 
as  well  as  other  data.  When  certain  simplifying  assumptions  are 
added  to  the  model,  the  equations  are  linearized  to  second  order 
in  deposited  energy,  and  the  deposition  profile  is  represented 
as  a  straight  line,  the  model  predicts  impulse  as  being  propor¬ 
tional  to  fluence,  multiplied  by  an  exponential  factor  that  con¬ 
tains  the  deposition  time.  That  is, 


I 


0.04186 


2Ck 


4>  exp 


where  I  is  the  impulse  in  ktap,  r  is  the  Gruneicen  ratio,  C.  is 

h  ^ 

the  bulk  sound  speed  in  cm/ysec  (K/p  )  ,  where  K  is  the  bulk 

2° 

modulus,  <J>  is  the  fluence  in  cal/cm  ,  r  is  the  electron  range  m 
cm,  and  tQ  is  the  deposition  time  in  ysec.  This  simplified  model 
gives  quite  good  agreement  with  magnesium,  aluminum,  and  silver, 
and  fair  agreement  with  lead,  using  no  adjustable  parameters  and 
handbook  values  for  material  properties. 


In  addition,  an  analysis  is  made  of  other  experimental 
techniques  for  measuring  impulse.  The  techniques  include  the 
use  of  quartz  gauges  for  measuring  stress  histories,  and  flyer 
plates  used  as  momentum  traps  on  the  back  of  tve  irradiated 
sample. 


Preceding  page  blank 
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B.  A  MODEL  FOR  MELT-DOMINATED  IMPULSE 

The  first  logical  step  in  building  a  model  for  melt- 
dominated  impulse  would  be  to  define  the  region.  A  quantitative 
estimate  of  the  upper  bound,  must  come  from  parameters  of  the 
model  itself.  For  now,  however,  the  assumption  will  be  made  that 
vapor  effects  are  negligible,  and  an  estimate  of  where  vapor 
effects  are  important  will  be  made  later  on. 


It  would  be  desirable  to  use  one  of  the  analytical  ex¬ 
pressions,  (Reference  17),  such  as  the  Whitener  (1)  or  BBAY  (2) 
expressions 

x 

I  =  \2  J  [EL  (m)  -  Ef  (m)]  ^  dm  (1) 

o 

x 

I  =  1.2  2  /  m  [EL  (m)  -  Ef(m)]  dm  %  (2) 

*o 

for  calculating  impulse.  To  do  so  requires  an  expression  for 
Ef (m) ,  the  final  internal  energy.  To  use  tne  sublimation  energy, 

E  ,  is  clearly  incorrect  since  this  would  imply  that  no  impulse 
is  generated  for  deposited  energy  less  than  Eg,  and  there  is  ex¬ 
perimental  evidence  that  impulse  is  generated. 


Certainly  the  McCloskey-Thompson  model  (Reference  17)  could 
be  used  to  calculate  the  final  state  energy.  When  this  is  done 
for  aluminum,  however,  the  model  does  not  fit  the  data.  Figure 
44  shows  the  impulse  data  in  aluminum,  and  the  impulse  predicted 
by  the  McCloskey-Thompson  model,  in  the  BBAY  expression,  where 
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2 

F]uence  (cal/cm  ) 


FIGURE  44,  ALUMINUM  IMPULSE  DATA  AND  McCLOSKEY-THOMPSON- 
BBAY  EXPRESSION 
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for  aluminum  E  i  -  160  cal/g  and  =  the  deposited  energy.  To 
understand  why  this  model  is  not  appropriate,  it  will  be  useful 
to  picture  the  process  going  on  in  melt-dominated  impulse  genera¬ 
tion  . 


Consider  energy  incident  from  the  left,  deposited  in  a 
material  instantaneously  and  for  simplicity,  with  a  linear  depo¬ 
sition  profile.  If  the  material  had  strength  and  were  a  linear 
fluid,  the  initial  stress  would  separate  into  halves,  a  right- 
moving  stress  wave  and  one  moving  to  the  left.  The  left-moving 
wave  reflects  from  the  front  surface  as  a  following,  anti-symmetric 
tension.  The  stress  history  at  a  point,  then,  shows  the  initial 
compression  (followed  by  an  increase  in  compression  if  the  point 
is  more  than  half  the  rarge  from  the  front  surface)  and  a  sudden 
drop  into  tension.  If,  on  the  other  hand,  the  heated  material 
has  no  tensile  strength,  an  increment  of  liquid  will  separate  from 
the  body  of  "solid*'  at  each  time  increment. 

In  effect,  the  vapor  phase  is  ignored  (the  equilibrium  vapor 
pressure  of  the  spalled  liquid  is  assumed  zero) ,  and  the  solid- 
liquid  transition  is  also  ignored — the  material  is  treated  as 
though  it  is  simply  a  strengthless  solid.  With  this  model,  the 
generated  impulse  can  be  calculated  by  finding  the  internal  energy 
changed  to  kinetic  energy  in  adiabatic  expansion,  and  this  energy 
substituted  in  the  Whitener  expression  for  impulse.  (The  Whitener 
rather  than  the  BBAY  expression  should  be  used  as  each  increment  of 
mass  is  suddenly  relieved  to  zero  pressure;  it  is  not  pushing  on 
other  mass.)  To  this  end,  the  kinetic  energy  for  a  simplified 
equation  of  state  is  derived  below. 

Assume  an  equation  of  state  to  be  of  the  form 

P  =  f(v)  +  ”  , 
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where  for  definiteness  the  units  of  pressure,  P,  and  internal 

12  3 

energy,  E,  are  10  erg/ref  cm  (Mbar) ,  v  is  the  relative  volume 
(V/V  -)  ,  and  I’  is  the  Gruneisen  ratio  (assumed  constant).  To 
find  the  path  of  adiabats,  the  condition  (dE/dv)  =  -P  is  substi¬ 
tuted,  to  give  the  differential  equation 


dP  +  r  +  1  p  _  1  d_ 

dv  v  v  dv 


which  has  the  solution 


P  =  f  (v)  -  r  (~)  f  C1  f  (C)  d£  +  C 

Specializing  f  to  be  the  first  term  of  a  standard 
equation  of  state  for  a  solid 

f<v>  -  S  (i  -  l) 

where  K  is  the  bulk  modulus. 


r  +  1 


polynomial 


-4) 


r  +  l 


K 


r  +  i 


and  C  is  determined  by  the  initial  condition 


rE 

po  ■  f  (vo>  +  77 


defining 


e  = 


P  = 


(r  *  i) 

K 


K 

r  +  i 


r  +  l 


(i  +  c)  i  ^1 


-  l 
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Thus  the  relative  volume  at  which  the  adiabat  is  at  zero  pressure, 
v*  is 


v*  =  v  (1  +  e) 


i/(r  +  l) 


The  energy  released  in  expanding  to  zero  pressure,  which  i~,  the 
kinetic  energy  of  motion,  is  given  by 


v’ 


AE  = 


/ 


v 


P(f-  <U 
adia 


Substituting 


AE  = 


v  (i  +  e)V(r  +  i) 

k  r° 

r  +  i  J 


V. 


o 


r  +  i 


(i  +  e) 


-  i 


dv 


Kv 


AE  = 


T(r  +  1) 


(f  +  1)  +  e  -  (T  +  1)  (1  +  e)1/(r  +  1) 


This  is  the  result  sought,  but  more  usefully  in  what  follows,  we 
expand  for  small  e  (using  a  Taylor's  series) ,  and  obtain 


Kv 


AE  = 


2  (r  +  i) 


v  p 

,2  ,  'o  Po 


2K 


This  expression  is  compared  to  the  McCloskey-Thompson  expression 
for  the  final  energy  in  Figure  45  for  the  parameters  of  aluminum 
assuming  the  energy  is  deposited  at  reference  volume.  It  is  clear 
that  the  final  energy  is  larger  for  the  linear  fluid;  hence,  the 
impulse  calculated  will  be  less  since  less  energy  is  available 
for  the  kinetic  energy  of  motion. 
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The  other  interesting  feature  of  the  expression  is  the 
dependence  on  the  volume  at  which  energy  is  deposited,  shown  in 
Figure  46,  for  the  parameters  of  aluminer  and  a  dose  of  0.1  Mbar 
900  cal/g) .  This  shows  qualitatively  that  if  the  deposition 
time  is  not  zero,  so  that  loading  near  the  front  surface  is  at 
relative  volumes  larger  than  one,  the  impulse  will  be  reduced. 

In  contrast  to  this,  the  expression  for  kinetic  energy  of  a 
Y -law  gas,  by  a  similar  derivation,  can  be  shown  to  be 

AE  =  <E0  -  Es> 

where  Es  is  the  unavailable  energy  (required  to  overcome  the 
potential  well  of  the  lattice) ,  independent  of  the  volume  at 
which  it  is  deposited.  Thus  we  have  the  first  indication  of  the 
importance  of  deposition  time  in  melt-dominated  impulse. 

With  the  approximations  of 

1.  Linear  deposition  profile 

2.  Material  removed  to  the  end  of  the  range,  r(cm) 

3.  Fluence  E  (1/2  c  p  r)  where  t  is  the  front  surface 

*o  Ko  o 

dose  in  (cal/g) 

4.  Linearized  expression  for  LE 

5.  Whitener  expression  for  impulse, 

the  impulse  is  given  by 


I  =  0.04186  * 

*  k 

where  I  is  the  impulse  in  ktap,  C,  is  the  bulk  sound  speed  in 

— - —  *  2 

cnt/usec  (=  /K/p  )  ,  and  <*>  is  the  fluence  in  cal/cm  .  This  is 

most  easily  derived  by  considering  that  the  momentum  density  (u) 

is  given  by 


This  assumes  that  bremsstrahlung  losses  and  electron  back- 
scattering  losses  are  neglected. 
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Inergy  Available  for  Motion  (10  erg/ref  cm 


1.0  1.1  1.2 


Relative  volume  during  instantaneous  deposition  (1) 
FIGURE  46.  ENERGY  AVAILABLE  FOR  MOTION  IN  ALUMINUM 
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where  AE  is  the  kinetic  energy  density  (10  erg/ref  cm  ) ,  u  is 
in  cm/usec,  and  integrating  the  velocity  over  the  mass  that  is 
included  in  the  electron  range. 

The  expression  for  impulse  is,  of  course,  only  appropriate 
for  nog]  I -|ible  deposition  times,  but  is  useful  for  predicting 
•ppor  bounds  for  melt-dominated  impulse.  Before  considering  the 
effects  of  deposition  time,  there  is  additional  information  con¬ 
tained  in  the  expression  for  impulse  and  kinetic  energy  that 
should  be  noted. 

First  the  expression  for  impulse  has  only  two  material 
sensitive  parameters,  T  and  C^,  and  for  metals,  the  ratio  varies 
by  about  a  factor  of  10  from  Be  (1.45)  to  Pb  (13.3).  Second,  the 
velocity  given  by  u  =  \J  (2AE/pq)  is  just  P/PQ  C^,  the  particle 
velocity  that  would  be  expected  from  the  Hugoniot  jump  conditions 
and  the  first  order  approximation  that  the  velocity  from  a  rare¬ 
faction  (the  Riemann  invariant)  Ur,  is  just  the  particle  velocity, 

U  .  Third,  since  neither  F  nor  C,  are  sensitive  to  small  alloy 
p  k 

concentration,  the  difference  in  impulse  between  different  alloys 
is  expected  to  be  small. 

Finally  the  linear  fluid  model  is  compared  with  the  McCloskey- 
->ompson  model  and  with  the  experimental  data  in  aluminum  (Figure 
47) 


The  effect  of  deposition  time  has  been  considered  by  Bade,  et 
al.,  (Reference  18),  in  their  calculation  of  non-blowoff  impulse 
1^.  When  their  formulas  are  generalized  to  an  arbitrary  deposition 
profile  (Reference  19),  the  expression  for  non-blowoff  impulse  is 


103 


c  t  r 

_L _0  f 

,2  2  J 

-k  CD  L  c  t  /2 


5k  V2 


c  1  (x) dx  + 


x  e  '  (x)dx 


When  this  formula  is  specialized  for  a  linear  deposition  profile, 
and  the  approximation  that  the  fluence,  4> ,  is  given  by 

*  =  k  eo  r  po 

is  retained,  the  expression  for  impulse  becomes 


I  =  0.0418  rfc-  <j>  f  (£) 
“ck 


-  C  +  1/3  ?  »  C  <  1 


where 


f(£)  = 


1/3  1/C 


e  2  1 


r  1  CkfcD 

c  -  2  — F~ 


CktD  5  3.5  r,  f(C)  s  e  f’  ±7% 

as  can  be  seen  in  Figure  48.  Thus,  to  the  accuracy  of  the  approxi¬ 
mations  we  may  use 

I  =  0.0418  ~  <b  e  -{CktD)/2r  ,  for  CktD  i  3.5  r 
k 


To  compare  the  expression  with  electron-beam  impulse  data  more 

readily.  Figure  2-13  of  Reference  20  has  been  used  to  obtain  the 

2 

range  of  electrons  in  g/cm  ,  pQr,  as  a  function  of  mean  electron 
energy  <&>  (MeV)  suitable  for  energies  near  0.25  MeV  for  all 


materials . 
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FIGURE  48 


COMPARISON  OF  DEPOSITION  TIME- DEPENDENCE 
FUNCTIONS 


From  this  figure 


Por  =  1.2  <g>2 


Hence  the  deposition  time  dependence  is  given  by 


f 


exp 


C.  t,.p 
k.  D*  o 


2.4 


'.€>2 


The  dependence  of  impulse  on  the  deposition  time  should  be 
strongly  emphasized.  For  melt-dominated  impulse,  it  becomes  an 
important  parameter.  If,  for  example,  an  attempt  were  made  to 
obtain  an  effective  E^  by  a  least-squares  fit  to  one  of  the 
analytical  expressions,  as  discussed  in  Reference  17,  and  the 
deposition  time  was  not  considered,  the  results  would  not  ex¬ 
trapolate  to  other  deposition  profiles  with  the  same  deposition 
time.  Indeed,  the  final  energy  Ef  could  be  a  strong  function  of 
deposition  time. 


Figures  4S,  50,  and  51,  taken  from  data  reported  in  References 
21  and  22,  show  impulse  versus  fluence  for  four  materials — Mg,  Al, 
Ag,  Pb — without  any  scaling,  scaled  by  the  linearized  impulse  ex¬ 
pression  without  regard  to  deposition  time,  arid  scaled  by  the 
full  formula,  respectively.  The  material  properties  used  are 
standard  PUFF  equation-of-state  parameters  and  are  shown  in 
Table  III.  From  these  figures,  when  the  approximate  time  de¬ 
pendence  and  impulse  formulas  are  used,  there  is  quite  reasonable 
agreement,  and  the  data  tend  to  cluster  about  the  theoretical  line 
that  has  no  adjustable  parameters. 

Since  the  expression  and  hence  the  model  for  melt-dominated 
impulse  has  been  verified,  we  may  use  the  model  to  estimate  the 
limits  of  the  melt-dominated  region. 
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TABLE  III 


PARAMETERS  USED  FOR  SOWING  IMPULSE  DATA 


r 

Ck 

(cm/usec) 

P°  3 
(g/cni  ) 

fcD 

(  VISOC) 

Al 

2.13 

0.536 

2.70 

0.03 

Mq 

1.46 

0.450 

1.74 

0.03 

Ag 

2.34 

0.327 

10.5 

0.03 

Pb 

2.78 

0.209 

11.3 

0.03 

Wc  assume  that  the  liquid  and  solid  both  have  the  same  bulk 

modulus,  K;  then,  the  intersection  of  the  release  adiabat  with  the 

liquid- vapor  dome  will  be  at  a  pressure  depending  on  the  initial 

internal  energy.  The  release  adiabat  of  the  liquid  that  intersects 

the  liquid  vapor  dome  at  about  0.1  kbar,  extended  to  zero  pressure, 

reaches  a  volume  of  approximately  1/2  v  ,  v  being  the  critical- 

c  c 

point  volume  (Reference  23).  The  0.1  kbar  level  is  chosen  more  or 
1  ess  arbitrarily  so  that  any  impulse  contribution  from  the  vapor, 
which  should  be  complete  by  85  4  usee,  is  less  than  0.2  ktap.  With 
the  approximation  that  vq  =  3vref»  we  estimate  the  internal  energy 
required  to  reach  this  adiabat  by  using 


v*  =  v  (1 
o 


+  e)1^  + 


If  the  energy  is  deposited  at  reference  volume,  this  gives 


/  %  _  Kx  10 

Eo  (cal/s)  -  4.186  p  r(r  +  1) 


#  -  -  f  +  1  . » 

(1.3  “  1) 


For  aluminum,  K  ==  0.75  Mbar,  T  =  2,  PQ  =  2.7  and  this  gives 


Eq  •»  2600  cal/g 

as  the  approximate  upper  bound  for  the  melt-dominated  impulse 
region . 
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C.  ANALYSIS  OF  OTHER  TECHNIQUES  USED  FOR  IMPULSE  MEASUREMENTS 

At  this  point, a  model  for  melt-dominated  impulse  has  been 
proposed,  and  a  simplified  formula  for  the  dominant  features  of 
impulse  generation  has  been  derived.  Tc  check  the  model  more 
carefully,  information  about  the  propagated  stress  wave  is  re¬ 
quired.  This  information  was  obtained  experimentally  by  measuring 
the  stress  history  with  quartz  gauges.  The  analytical  details  of 
stress  propagation  were  obtained  by  programming  the  model  for  liquid 
splashoff  as  part  of  the  constitutive  relation  in  a  Lagrange  hydro- 
dynamic  code,  POD. 

The  solid -liquid  phase  transition  was  ignored.  The  justification 
for  this  is  twofold:  (1)  The  results  of  McQueen,  et  al.  (Reference  24) 
show  that  the  liquid-solid  transition  in  copper  was  virtually  indis¬ 
tinguishable  on  a  shock  speed-particle  velocity  plot,  and  (2)  that 
inclusion  of  the  phase  transition  was  unnecessary  to  obtain  good 
agreement  with  stress  and  impulse  data.  The  constitutive  relation, 
then,  describes  the  liquid  phase  as  having  the  same  bulk  properties 
as  the  solid,  but  with  yield  stress,  shear  modulus,  and  tensile 
strength  that  approach  zero  as  the  internal  energy  reaches  the  en¬ 
thalpy  required  to  reach  the  melting  point  at  1  atm.  The  relations 
used  are  shown  in  Table  IV  and  Figure  52.  With  these  constitutive 
relations,  machine  diagnostics  were  used,  as  described  in  Reference  2, 

to  obtain  the  energy  spectrum,  mean  angle  of  incidence  of  the  electron 

stream,  fluence,  and  effective  deposition  time.  The  deposition 
profile  was  then  calculated  by  a  Monte  Carlo  electron  transport  code, 

c*  i  nn  wac*  i'  r,  p  c<ar  ^  oc?  r»f  V\\r  rJ  prri  i  r*  r*nHp 

V*  A  ft  >-*  Wtl  X  X  •  •  V  <>  MIM  X  V  ft  ft  *»  V*  S>  M  U  \>  Vft  X  ft  ft  ^4  ft)  WX  V*  O  O  X  ftftj  ftftx  W  J  ft  ft  Mill  X  *«• 

calculations . 
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TABLE  IV 


CONSTITUTIVE  RELATION  FOR  ALUMINUM 

Pressure 

P  =  0 . 75y  +  0 . 4 0 5 )i2  +  (2  +  2p)  E  ,  u  ^  1  (Ref  26) 

0.75p  +  (2  +  2p)  E  ,  u  <  1 

Where  u  .1  p/po~l 
Shear  Modulus 

G  =  0.275  -  15.2SE  ,  E  <  0.0181  (Ref  26) 

0  E  >  0.0181 

Von  Mises  yield  stress,  YQ  =  f{£),  f (E)  is  shown  in  Figure  52. 
Spall  strength  P  .  (E)  =  -0.006  (1  -  E/0.0181). 
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The  first  series  of  calculations  was  made  to  obtain  impulse 
by  recording  the  momentum  of  the  spalled  liquid.  The  results  of 
the  calculations  and  the  simplified  model  discussed  in  Section  III 
are  shown  with  the  impulse  data  in  Figure  53.  In  addition  to  these 
calculations,  other  calculations  were  made  to  verify  the  dependence 
of  impulse  on  various  constitutive  parameters.  They  showed  that, 
indeed,  a  20%  reduction  in  bulk  modulus  (programmed  as  a  linear 
decrease  with  energy)  caused  a  10%  increase  in  impulse,  and  that 
the  effect  of  deposition  time  could  be  represented  by  exp  (-  t^/ 

2r)  . 


The  second  series  was  run  to  obtain  the  predicted  stress 
histories  in  quartz.  The  deposition  profile  calculated  for  one  of 
these  is  shown  in  Figure  54 .  To  obtain  agreement  with  the  measure¬ 
ment,  it  was  found  necessary  to  include  the  epoxy  bond  between  the 
quartz  and  the  aluminum  in  the  calculations.  When  this  was  done, 
the  experimental  and  calculated  stress  histories  were  compared  and 
are  shown  in  Figure  55 .  The  good  agreement  lends  considerable 
credibility  to  the  proposed  model. 

When  the  quartz  gauge  record  is  studied  more  carefully,  it  is 
noted  that  the  integral  of  stress  over  time  is  1.3  ktap,  which  is 
less  than  the  calculated  momentum  in  the  spalled  liquid,  which  was 
1.6  ktap.  The  impedance  mismatch  of  quartz  and  aluminum  indicates 
that  if  the  stress-time  momentum  in  alur inum  were  1.6  ktap,  the 
quartz  would  show  at  least  1.5  ktap.  Thus,  we  have  a  small  but 
apparently  real  discrepancy  in  the  calculation.  To  understand  this 
discrepancy,  one  must  consider  the  details  of  elastic-plastic 
behavxor . 
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Impi 


FIGURE  53.  ELASTIC-PLASTIC  HYDRODYNAMIC 
CALCULATIONS  OF  IMPULSE 


Calculation  With  1  mil  Epoxy  Bond 
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FIGURE  55.  STRESS  HISTORY  IN  REAR-SURFACE  QUART/.  GAUGE 


If  enough  material  at  the  front  surface  is  removed,  the  normal 
tensile  tail  caused  by  thermal  stress  production  (Reference  27)  is 
completely  suppressed,  and  the  traveling  stress  pulse  is  only  com¬ 
pressive.  'Ihe  path  in  stress-volume  space  that  is  traversed  by 
elastic-plastic  compression  and  release  has  hysteresis;  that  is, 
there  is  a  residual  compression.  This  same  hysteresis  occurs  in 
stress-particle  velocity  space,  so  that  the  result  of  passing  a 
co  tiprossive  stress  wave  is  to  leave  the  material  at  a  uniform 

velocity.  This  velocity,  for  the  parameters  of  aluminum  used, 

-4 

corresponds  to  5  x  10  cm/Msec,  which  is  0.135  ktap/mm.  The 
distance  that  the  stress  wave  traveled  to  reach  the  guartz  in  the 
above  calculation  was  greater  than  1  mm,  which  accounts  for  the 
discrepancy.  This  phenomenon  will  occur  with  ail  materials  that 
have  different  loading  and  unloading  paths.  Further,  if  the  plastic 
part  of  the  stress  wave  in  the  irradiated  sample  is  transmitted  as 
a  stress  wave  that  stays  within  the  elastic  limit  of  quartz,  the 
residual  velocity  of  the  material  will  appear  as  a  stress  offset 
in  the  quartz,  an  apparent  baseline  shift. 

A  second  technique  suggested  for  the  measurement  of  impulse  is 
the  use  of  flier  plates  attached  to  the  rear  of  irradiated  samples. 
The  velocity  of  these  plates  can  be  measured  by  recording  their 
trajectories  with  a  high-speed  movie  camera,  and  measuring  the  mass 
and  center-of-mass  velocity.  To  .relate  this  momentum  to  the  ejecta 
momentum,  we  will  consider  the  one-dimensional  stress-wave  inter¬ 
actions  in  some  detail. 

Assume  that  the  flier  and  target  have  the  same  shock  impedance, 
and  that  the  stress  wave,  however  generated  in  the  target,  is  an 
elastic  finite  pulse.  The  first  important  feature  is  that  the  simple 
energy  and  momentum  conservation  relations  used  for  rigid  bodies  do 
not  work.  That  is,  in  the  lab  frame, 


1,  > 


when  the  thickness  of  the  pi.* to,  F,  is  such  that  the  entire  stress 
pulse  can  be  contained  in  it. 

This  flier  velocity  is  certainly  a  contradiction,  and  at  first 

thought,  it  appears  to  be  nonconservation  of  kinetic  energy.  That  is, 

2  2 

mFVF  ^  mTVT  *  The  kinetic  energy  that  is  apparently  lost, 

however,  gees  into  the  ringing  of  the  elastic  pulse,  which  becomes  ap¬ 
parent  when  it  is  remembered  that  for  a  nonuniform  velocity  distribution — 
a  pulse — the  square  of  the  average  velocity  is  less  than  the  average 
of  the  squared  velocity.  That  is.  when  the  kinetic  energy  associated 
with  the  momentum  of  a  stress  pulse  is  assumed  smoothed  into  conter- 
of-mass  translation,  kinetic  energy  is  lost. 

When  the  flier  has  higher  shock  impedance  than  the  target,  and 
the  flier  is  thick  enough  to  contain  the  pulse  for  elastic  stresses, 


2Zp  (mTvT) 

VF  “  Z"T  +  Zp  in^ 

where  Zc,  and  Z,„  are  the  respective  shock  impedances.  If  the  impedance 

r  X 

of  the  flier  is  less  than  the  sample 


(^tvt)  /  2Zp  \ 
mp  y  T  ^F  ) 
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how  much  less  v„  is  than  the  upper  bound  depends  on  the  shape  of 
the  stress  pulse. 


If  there  is  a  hysteresis  in  the  sample,  then  the  flier,  like 
the  quartz  gauge,  will  trap  only  part  of  the  momentum.  For  "thick” 
fliers,  with  shock  impedances  greater  than  or  equal  to  the  target 


21, 


mFVF 


(iVt  “  po,ur) 


whe--*e  ur  is  the  residual  sample  velocity.  (For  aluminum  the  sub¬ 
tractive  term  is  *=  0.13  ktap/mm.) 


'i’hings  are  not  nearly  so  simple,  however,  when  the  flier  is 
not  thick  enough  to  contain  the  stress  pulse,  due  to,  say,  vapor 
pushing  for  long  tiroes.  In  this  case,  a  calculation  of  the  stress 
waves  is  ’required  to  get  even  an  approximate  value  for  the  momentum 
fraction  trapped  in  the  *lier. 

D.  EJECTA  VELOCITIES 

A  further  test  of  the  spall-dominated  model,  and  probably  more 
direct  confirmation  than  even  transmitted  stress,  would  be  agreement 
with  calculated  and  measured  ejecta  velocities. 

Under  a  program  carried  out  at.  Physics  International  by  D,  Dean 

of  the  Sandra  Corporation,  (Reference  9},  a  number  of  experiments  on  the 

Pulserad  733  machine  were  performed  using  a  high-speed  framing  camera  to 

record  the  motion  of  the  front  surface  ejecta.  Selected  frames  of 

2 

the  front  surface  response  of  aluminum  exposed  to  40  cal/cro  of  0.26 
MeV  electrons  are  shown  in  Figure  36.  The  position  leading  and  trail¬ 
ing  edges  of  persistent  structure  in  the  ejecta  show  velocities  of 
0,06  and  0.020  cm/tsec.  The  calculated  mean  velocity  for  the  deposi¬ 
tion  conditions  was  0,03  cm/ysec .  It  is  seen  that  good  quantitative 
agreement  between  calculated  and  measured  velocities  exists  within 
the  uncertainty  of  the  experiment,  thereby  providing  confirmation  ol 
the  melt-dominated  model. 
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SECTION  V 


CONCLUSIONS 


Measurements  of  impulse  generation  using  a  ballistic  pendulum, 
transmitted  stress  histories  using  quartz  gauges,  and  ejecta  ve¬ 
locities  using  framing  cameras  have  been  successfully  performed  in 
an  electron-beam  environment.  These  measurements,  together  with 
analysis,  have  been  used  in  an  extensive  program  to  study  impulse 
generation  in  aluminum  in  the  melt-dominated  regime. 


The  major  result  of  the  combined  analytical  and  experimental, 
program  is  the  conclusion  that  the  description  of  material  removal 
by  liquid  spall  extends  to  rather  large  energy  densities.  In 
aluminum,  for  example,  it  is  anticipated  that  energy  densities 
nearly  sufficient  to  cause  complete  vaporization  at  atmospheric 
pressure  can  be  successfully  treated  by  ignoring  the  vapor  phase. 
This  analysis  was  also  seen  to  be  a  reasonable  approximation  for 
other  simple  materials.  Thus,  we  take  it  as  a  general  conclusion 
that  the  liquid-spal]  model  is  appropriate  up  to  energy  densities 
given  approximately  by 


Eq  (cal/g) 


K  x  10' 


4.186  pQ  nr  +  1) 


(1.5 


r  + 


-  1) 


where  K  is  the  bulk  modulus  in  Mbar,  T  is  the  Gruneisen  ratio,  and 

3 

PQ  is  the  initial  density  in  g/cm  .  When  this  formula  is  evaluated 
for  various  materials.  Mg,  Al,  Ag,  and  Pb,  the  maximum  dose  is 
within  25%  of  the  energy  of  vaporized  metal  at  atmospheric  pressure. 


With  recent  improvements  in  achievable  doses  on  the  Pulserad 
728,  it  is  anticipated  tha*.  an  extension  of  the  experimental  results 
presented  here  would  be  made  to  experimentally  verify  the  upper 
bound  of  the  spall-dominated  impulse  region  and  to  determine  appro¬ 
priate  models  in  the  mixed  phase  and  vapor-dominated  regions.  A 
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combined  analytical  and  experimental  approach,  such  as  used  on  this 
program  and  outlined  in  Figure  57,  is  recommended.  In  this  way,  the 
results  ol  experiments  and  analysis  can  affect  each  other  during 
the  course  of  the  program. 

It  is  expected  that  such  an  approach  will  provide  the  funda¬ 
mental  basis  of  understanding  required  to  model  mixed-phase  effects 
and  permit  the  extension  to  even  more  complex  materials  such  as 
fiber-reinforced  composites. 
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